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Preface 


The Gamma function is an extension of the factorial func- 
tion, with its argument shifted down by 1, to real and 
complex numbers. The Gamma function is defined by 
an improper integral that converges for all real numbers 
except the non-positive integers, and converges for all 
complex numbers with nonzero imaginary part. The fac- 
torial is extended by analytic continuation to all complex 
numbers except the non-positive integers (where the inte- 
gral function has simple poles), yielding the meromorphic 
function we know as the Gamma function. 


The Gamma function has very many extremely impor- 
tant applications in probability theory, combinatorics, 
statistical and quantum mechanics, solid-state physics, 
plasma physics, nuclear physics, and in the decades-long 


quest to unify quantum mechanics with the theory of 
relativity — the development of the theory of quantum 
gravity — the objective of string theory. 


The problem of extending the factorial to non-integer ar- 
guments was apparently first considered by Daniel Bernoulli 
and Christian Goldbach in the 1720s, and was solved at 
the end of the same decade by Leonard Euler. 
Euler gave two different definitions: the first was an infi- 
nite product, of which he informed Goldbach in a letter 
dated October 13, 1729. He wrote to Goldbach again on 
January 8, 1730, to announce his discovery of the inte- 
gral representation. Euler further discovered some of the 
Gamma function’s important functional properties, no- 
tably the reflection formula. 


Carl Friedrich Gauss rewrote Euler’s product and then 
used his formula to discover new properties of the Gamma 
function. Although Euler was a pioneer in the theory of 
complex variables, he does not appear to have consid- 
ered the factorial of a complex number, as Gauss first 
did. Gauss also proved the multiplication theorem of 
the Gamma function and investigated the connection be- 
tween the Gamma function and elliptic integrals. 


Karl Weierstrass further established the role of the Gamma 
function in complex analysis, starting from yet another 
product representation. Weierstrass originally wrote his 
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product as one for 1/T, in which case it is taken over 
the function’s zeros rather than its poles. Inspired by 
this result, he proved what is known as the Weierstrass 
factorization theorem — that any entire function can be 
written as a product over its zeros in the complex plane; 
a generalization of the fundamental theorem of algebra. 


The name of the Gamma function and its symbol I’ were 
introduced by Adrien-Marie Legendre around 1811; Leg- 
endre also rewrote Euler’s integral definition in its mod- 
ern form. The alternative "Pi function" notation H(z) = 
z! due to Gauss is sometimes encountered in older lit- 
erature, but Legendre’s notation is dominant in modern 
works. It is justified to ask why we distinguish between 
the "ordinary factorial" and the Gamma function by us- 
ing distinct symbols, and particularly why the Gamma 
function should be normalized to T(n + 1) = n! instead 
of simply using "I (n) = n!". Legendre’s motivation for 
the normalization does not appear to be known, and has 
been criticized as cumbersome by some (the 20th-century 
mathematician Cornelius Lanczos, for example, called it 
"void of any rationality" and would instead use z!). Leg- 
endre’s normalization does simplify a few formulas, but 
complicates most others. 


A large number of definitions have been given for the 
Gamma function. Although they describe the same func- 
tion, it is not entirely straightforward to prove their equiv- 
alence. Instead of having to find a specialized proof for 


each formula, it would be highly desirable to have a gen- 
eral method of identifying the Gamma function given any 
particular form. 


One way to prove equivalence would be to find a differ- 
ential equation that characterizes the Gamma function. 
Most special functions in applied mathematics arise as 
solutions to differential equations, whose solutions are 
unique. However, the Gamma function does not appear 
to satisfy any simple differential equation. Otto Holder 
proved in 1887 that the Gamma function at least does 
not satisfy any algebraic differential equation by show- 
ing that a solution to such an equation could not satisfy 
the Gamma function’s recurrence formula. This result is 
known as Holder’s theorem. 


A definite and generally applicable characterization of 
the Gamma function was not given until 1922. Harald 
Bohr and Johannes Mollerup then proved what is known 
as the Bohr-Mollerup theorem: that the Gamma function 
is the unique solution to the factorial recurrence relation 
that is positive and logarithmically convex for positive z 
and whose value at 1 is 1 (a function is logarithmically 
convex if its logarithm is convex). 


The Bohr-Mollerup theorem is useful because it is rela- 
tively easy to prove logarithmic convexity for any of the 
different formulas used to define the Gamma function. 
Taking things further, instead of defining the Gamma 
function by any particular formula, we can choose the 


conditions of the Bohr-Mollerup theorem as the defini- 
tion, and then pick any formula we like that satisfies the 
conditions as a starting point for studying the Gamma 
function. This approach was used by the Bourbaki group. 


G.P. Michon describes the Gamma function as "Arguably, 
the most common special function, or the least 'special! 
of them. The other transcendental functions ... are called 
'special' because you could conceivably avoid some of 
them by staying away from many specialized mathemat- 
ical topics. On the other hand, the Gamma function is 
most difficult to avoid." 


The Gamma function finds application in diverse areas 
such as quantum physics, statistical mechanics and fluid 
dynamics. The Gamma distribution, which is formulated 
in terms of the Gamma function, is used in statistics to 
model a wide range of processes; for example, the time 
between occurrences of time-series events. The primary 
reason for the Gamma function’s usefulness is the preva- 
lence of expressions of the type f(t) exp(—g(t)) which 
describe processes that decay exponentially in time or 
space.’ Integrals of such expressions can often be solved 
in terms of the Gamma function when no elementary so- 
lution exists. For example, if f is a power function and 
g is a linear function, a simple change of variables yields 


+ See chapter 2 for further discussion of this topic. 
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The fact that the integration is performed along the en- 
tire positive real line might signify that the Gamma func- 
tion describes the cumulation of a time-dependent. pro- 
cess that continues indefinitely, or the value might be 
the total of a distribution in an infinite space. It is 
of course frequently useful to take limits of integration 
other than 0 and oo to describe the cumulation of a fi- 
nite process, in which case the ordinary Gamma func- 
tion is no longer a solution; the solution is then called 
an incomplete Gamma function. (The ordinary Gamma 
function, obtained by integrating across the entire posi- 
tive real line, is sometimes called the complete Gamma 
function for contrast.) 


The Gamma function’s ability to generalize factorial prod- 
ucts immediately leads to applications in many areas of 
mathematics; in combinatorics, and by extension in areas 
such as probability theory and the calculation of power 
series. Many expressions involving products of successive 
integers can be written as some combination of factorials, 
the most important example perhaps being that of the 
binomial coefficient. 


By taking limits, certain rational products with infinitely 
many factors can be evaluated in terms of the Gamma 
function as well. Due to the Weierstrass factorization 
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theorem, analytic functions can be written as infinite 
products, and these can sometimes be represented as fi- 
nite products or quotients of the Gamma function. For 
one example, the reflection formula essentially represents 
the sine function as the product of two Gamma functions. 
Starting from this formula, the exponential function as 
well as all the trigonometric and and hyperbolic functions 
can be expressed in terms of the Gamma function. 


The hypergeometric function and special cases thereof, 
can be represented by means of complex contour inte- 
grals of products and quotients of the Gamma function, 
called Mellin-Barnes integrals. 


The Gamma function can also be used to calculate the 
"volume" and "area" of n-dimensional hyperspheres. 


An elegant and deep application of the Gamma func- 
tion is in the study of the Riemann zeta function. A 
fundamental property of the Riemann zeta function is 
its functional equation. Among other things, it provides 
an explicit form for the analytic continuation of the zeta 
function to a meromorphic function in the complex plane 
and leads to an immediate proof that the zeta function 
has infinitely many so-called "trivial" zeros on the real 
line. Borwein et. al call this formula "one of the most 
beautiful findings in mathematics". 
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The Gamma function has caught the interest of some 
of the most prominent mathematicians of all time. In 
the words of Philip J. Davis, "each generation has found 
something of interest to say about the Gamma function. 
Perhaps the next generation will also." Its history re- 
flects many of the major developments within mathemat- 
ics since the 18th century. 
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Analytic Continuation of the 
Factorials 


Analytic continuation is a technique used to extend the 
domain of a given analytic function. Analytic continua- 
tion often succeeds in defining further values of a func- 
tion. Consider the sum of the first n natural numbers: 
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We may simply construct a new formula, not involving 
summation, in the following way: 


1+2+---+(n-1)+n 
n+(n—1)+---+24+1 
(n+1) +--+ (n+1) 


n terms 


Gauss famously discovered this independently in his early 
childhood. In doing this, we have performed the sum 
twice, so 


s= n(n Fi (1.1) 


It is not meaningful to speak of the first 3/2 natural 
numbers, for instance, but the above equation does in- 
terpolate between the integer values. We may apply the 
function to non-integer values and get, for example, 

3 1 3 5 15 

aenn S, 

2 2 2 2 8 

The factorial function n! is also only defined over the 
positive integers: 


nt :=n-(n—1)-(n-2)----- 3-2-1 (1.2) 


with the convention that 0! = 1. A formula such as 
Equation 1.1 which would allow computation of n! with- 
out having to perform the multiplication was sought by 
the mathematics community for many years. 


Euler first discovered a formula for the factorials that 
allows us to compute a factorial without performing the 
multiplication — the so-called elementary definition of the 
Gamma function, is Euler’s integral, which is as follows: 


Definition 1.1. (Gamma function). For 0 < z < o, 
xeER, 


r(x) := [ Gore dt. (1.3) 


The defining improper integral itself, in the real case, con- 
verges if and only if 0 < x < co. This integral function 
is extended by analytic continuation to all real numbers 
except the non-positive integers, x ¢ {{0} UZ}, by the 
use of a recurrence formula to yield the Gamma function. 
The question may be asked as to how the Gamma func- 
tion is computed on the negative nonintegers given that 
the defining integral does not converge on them. Indeed, 
Euler’s integral does not converge for x < 0, but the 
function it defines over the positive reals has a unique 
analytic continuation to the negative reals. One way to 
find that analytic continuation is to use Euler’s integral 
for positive arguments and extend the domain to negative 
numbers through the application of the identity, 


T(z +n) 


Deas (ay 


T(x) 


4 -2 2 4 


Figure 1.1: The Gamma function plotted over a portion of the real 
line. 


choosing n sufficiently large such that z + n is positive. 
The product in the denominator is zero when z equals any 
of the integers 0, —1, —2,.... Thus, the gamma function 
must be undefined at those points. 


A restriction on the domain of the Gamma function does 
not equate to a restriction on the range of integration 
over the variable of integration (the argument t of the in- 
tegrand). But rather, the argument of the Gamma func- 
tion parametrizes the integrand, say g,(t) := t®~te~‘, 
given in the Gamma function’s definition — restrictions 
on the domain of the Gamma function restrict the al- 
lowed values of the parameter x defining g,(t). 
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Figure 1.2: The argument of the Gamma function parametrizes the 
curve representing the integrand gz(t) in the Gamma function defini- 
tion. It is a single curve, such as one of the curves shown here, which is 
integrated over the range of integration (from t = 0 to oo) in the eval- 
uation of the Gamma function. Here are shown plots of the integrand 
gx(t) = e~'t®~! for parametrizations x = 1,2,3 and 4. If one were to 
smoothly change the value of the parameter x from 1 to 4 then the shape 
of the curve would smoothly change from that of the x = 1 curve to that 
of the x = 4 curve. 


We shall prove a result known as the functional equation 
(a recurrence relation) later in this chapter. It is defined 
as 


T(z+1) = 2I(z). 


Inverting the functional equation we have 
1 
T(z) = -P(z4+1). 
Z 
When z = 0 it diverges because 


r(1) =. edt = lim [-e“]> = lim (~e™+1) =1 
0 


m—> co m—-> oo 


is finite. 


Formally, we can discover that the Gamma function has 
simple poles in all negative integers by simply iterating 
the recursion: 


Peay = E, 


2, Ae r0) (n EN). 
n! 

The power function factor that appears within the inte- 
grand of the Gamma function definition may have a com- 
plex exponent. What does it mean to raise a real number 
to a complex power? Ift is a positive real number, and z 
is any complex number, then t?~! is defined as e@~) 8 , 
where x = log(t) is the unique real solution to the equa- 
tion e” = t. By Euler’s formula, e* = cos@ + isin 0, so 
we have for example: 


g3+4i L 93,94 _ gefilog2 
= 8(cos(log 2) + i- sin(log 2))* 
—7.4615 + 2.88549: 
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Though we can straightforwardly compute the complex 
power of a real number in the analytical sense, a more 
intuitive interpretation of taking the complex power of a 
real number is elusive. Nonetheless, by virtue of the re- 
currence relation, T (z) can be interpreted as a generaliza- 
tion of the factorial function to complex numbers. Since 
T(z) =T(z+1)/z and I'(z+1) is defined for Re(z) > —1 
we can extend analytically the domain of T (z) to the strip 
—1 < Re(z) < 0. Repeating the same argument, we can 
extend I (z) to all values of the complex plane, except for 
B= Ot a1 EE 


Since the argument of the Gamma function may be com- 
plex, we may restate the definition of the Gamma func- 
tion in the following way: 


Definition 1.2. (Gamma function). For |z| < o, 
z €C, 


E= [ te “de. (1.4) 


The improper integral itself converges iff Im(z) 4 0, or 
z€Rand0<z<o. The integral function is extended 
by analytic continuation to all complex numbers except 
the non-positive integers, z € {{0} UZ}, through the 
use of a recurrence formula yielding the Gamma function. 


Let’s take a closer look at definition 1.2. Notice that we 
are integrating over the range [0, 00) — here we have an 


Figure 1.3: A surface plot of the absolute value of the Gamma function 
over a portion of the complex plane that illustrates the pattern of its 
domain of convergence, defined as z € C such that z € R and z ¢ {{O}U 
Z`}, or Im(z) £0. 


example of an improper integral, which are defined by a 
limiting process: 
i Hoars a] Fae 
noo m 

If Im(z) = 0, then when z < 1, we have that is 
undefined at t = 0, so we must, in that case, define our 
improper integral in terms of two limits: 

n—oo m—-0 


I tle dt := lim lim ite tdt 
0 


m 


Next, let’s verify that the integral in definition 1.2 con- 
verges. To this effect, we shall need the following lemma: 


Lemma 1.1. Comparison Test for Improper In- 
tegrals. If f(x) > g(x) > 0 on the interval [a, oo) then, 


o If [7 f(x) dx converges, then so does f? g(x) da. 
o If [7 g(x) dx diverges, then so does [~™ f(x) da. 


We shall take this as intuitively obvious and not pursue 
a formal proof here. 


Theorem 1.1. f t?~'e~' dt converges for all Re(z) € 
(0, 00). 

Proof. Let x = Re(z) and split the integral into a sum of 
two terms: 


oo 1 oo 
| tle dt = I tte dt + | tle dt. 
0 0 1 


The strategy is to show that each of the two integrals 
on the right-hand side converge — in both cases we con- 
struct f(t) such that 0 < ¢”-le* < f(t) and show that 
the integral of f(t) over (0,00) is finite, and hence so is 
the integral of t?~te~*. 


In the case of the first integral, since e™% < 1 for t > 0, 
we have that 0 < t®~te~* < t®-! for t > 0. Then 


1 1 
o< f metats | t! dt 
0 0 


ale 1 a 
= lim |— = lim |--— |. 
a>0t+ |£ |a GOR z 
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If x > 0, then aë — 0 as a — 0+, so that b t1 dt 
converges to +. Hence We t?~1e—t dt converges. 


In the case of the second integral, first note that ttet > 
0 as t > œ for any q € R. Hence for any x € R, there 
exists € such that 0 < t*~!e~*/2 < 1 fort = e. So we 
further split the second integral at e: 


| pid= metas f ite dt. 
1 1 € 


The first term is finite, being a finite integral. For the 
second term, when t > € we have 


thus 


/ ite" dt < | e™ dt = lim [-2e7#/?]" 


a—oo t=e 


= lim (2e~*/? — 2e-°/?) = 267°, 
Q— OO 
This shows the second term is convergent, thus [P7 t?~'e~ dt 
converges. All said, Wee t?-1e— dt converges for all x such 
that 0 < x < oo by lemma 1.1. 


Next, it is proper to ask what the definition of T(x) 
means. It has already been said that T(x) extends the 
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factorial function to 0 < x < oo, but we have not shown 
this — in order to do so, we need the following result: 


Proposition 1.2. For all x > 0, T(x + 1) = r(x). 


Proof. From definition 1.1, 
T(i¢+1)= I Pe "de: 
0 

Integration by parts’, with f(t) = t® and g(t) = —e™, 
yields 

f te dt — f atle dt = [-t*e"']>. 

0 0 
Rearranging, 


Tr(x+1)= lim [wey + f rt 26 dt 
0 


m—->oo 0 


= lim (~e m”) +a f le” dt. 
m00 0 


+ The product rule says that 
(f-g =f- 9 +F g, 


from which it follows that 


b b 
b 
f toas=- f f-o acto lt, 
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Now, as m — co, m*e—™ — 0, so we have 


T(x +1) = of (te. “dt = rT (x). 
0 


The recurrence relation, also known as the functional 
equation, 


I(z+1) = 2IT(z) 


is valid for all complex arguments z for which the integral 
in definition 1.2 converges. 


We may extend the functional equation to negative values 
through inversion, 


r 1 
T(z) = ad), 
Z 
so for example [(—1/2) = —2T(1/2). Reiteration of 


this identity allows us to define the Gamma function on 
the whole real axis except on the nonpositive integers 
(0i = N 


That said, we may now show the following: 


Theorem 1.3. For alln EN, T T(n+1)=n!. 
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Proof. From proposition 1.2, we have 


Tin +1) =n-T(n) 
=n-(n—1)-T(n-1) 


It remains to show that ['(1) = 1 = 0!. By definition 1.1, 
we have 


(1) =i edt = lim [-e“]" = lim (—e™+1) = 1. 
0 


mM—> o0 m— o0 


Euler wrote a letter on January 8, 1730 to Christian 
Goldbach in which he proposed the following definition 
for the Gamma function: 


Definition 1.3. Let x > 0, 


T(x) = J (—log(t))*-! dt (1.5) 


We can show this historical definition equivalent to its 
more usual form, definition 1.1, through a simple change 
of variable u = — log(t). 
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Proof. Let u = —log(t). Then we have 


dt 


— and du = ——. 


t=e 
Our corresponding limits of integration are 
t=0—>u=œ and t=1—>u=0. 


We thus have 


0 oo 
-f u” te™ du = f u® te“ du 
oo 0 


where a simple relabeling yields definition 1.1. 


Other forms of the Gamma function are obtained through 
a simple change of variables, as follows: 


r(z) = 2 | ye te¥ dy by letting t = y? 
0 


1 1 z—1 
r(z) = n (in z) dy by letting e™% = y. 
0 


Definition 1.4. Derivatives of the Gamma Func- 
tion. The derivatives of the Gamma function can be 
calculated by straightforward differentiation under the 
integral sign: 


I’ (x) = tte‘ log(t) dt 


0 


r(x) = [ tle log" (t) dt 
0 
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where the identity 4a” = a” log(a) is used’. 


So how do we actually compute the Gamma, function? 
It was previously mentioned that we can compute the 
Gamma function for negative argument using the recur- 
sion formula: 


[(z+n) 


PEIS z(z+ 1) (z+n+1) 


where we choose n sufficiently large such that z + n is 
positive. Another formula for computing the Gamma 
function for negative argument involves Euler’s reflection 
formula (discussed later): 


zT (z)sinrz 


However both of the above two formulas rely upon an 
ability to compute the Gamma function for positive ar- 
gument. 


The integral definition is not very useful in terms of effi- 
ciency; to produce an accurate result, an extremely high 


f This derivative identity follows from a simple application of 
the chain rule: 
d x£ d zloga _ ploga* 


T’ A -loga =a" - loga. 
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number of terms would have to be added during some 
numerical integration procedure. There are several other 
ways, possessing various degrees of efficiency, to numeri- 
cally compute the Gamma function. The simplest to un- 
derstand is as follows: A series expansion (derived later) 
for log(T (1 + x)) exists such that for |x| < 1 we have 


log(T(1 + £)) =-y2+ 5° ca (1.6) 


It follows easily from the equation above and the func- 
tional equation that 


tl (=1)*¢(k)z" 
coy: z exp (x-5 ez) (1.7) 


k=2 


Other formulas for computing the Gamma function in- 
clude the Lanczos approximation, which is used to com- 
pute the Gamma function to arbitrary precision, and 
Stirling’s formula, which is an asymptotic formula used 
to approximate the value of the Gamma function given 
very large arguments when absolute precision is less im- 
portant, and is discussed in a later chapter. 


Integrals Involving a 
Decaying Exponential 


We begin with the Gamma function, which is defined by 
Euler’s integral: 


T(z) =| dzz” te, for Re(z) > 0. 
0 


It can be immediately recognized that if we set z = n+ 1 
then we have 


Tin+1)= f tet dt, for Re(n) > —1. 
0 


We arrive at a useful relation if we make the substitution 
x = vt such that dx = vdt in the evaluation of the 
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following integral: 


I dif'e 
0 


de £N” oe 
[ FQ" 
fe gre "dz T(n+1) 


prtl n prtl 


Integrals over temporally or spatially decaying processes 
(such as collisional damping at rate v ~ 1/7) often result 
in integrals of the form 


l dtt?e/7 = oy, dxx"e~” where x = t/r. 
0 0 


Two values of the argument z of fundamental interest 
for Gamma functions are z = 1 and z = 1/2. For z=1 
the Gamma function becomes simply the integral of a 
decaying exponential: 


r(1) = f dre” =1. 
0 


For z = 1/2, by using the substitution z = u? the Gamma 
function becomes the integral of a Gaussian distribution 
over an infinite domain: 


1 oa 2 X 2 
Tr(+)= 2f du ufu?) e7 = 2f due“ = yT. 
0 0 


+ Refer to Chapter 5 for a proof of this result. 
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When the argument of the Gamma function is a positive 
integer (given by the map z > n > 0), the Gamma 
function simplifies to a factorial function: 


T(in4+1) =nI(n) =---=n(n—1)(n— 2)---L= nl. 


Using this factorial form for the Gamma function, one 
thus finds that 


f dtte ™ = r°ttn!, forn = 0,1,2,... 
0 


using the usual convention that 0! = 1. The first few of 
these integrals are 


1 1 1 
f — t/r pe = | dr4 ape a we 1 
0 T | Er 0 x? 2 


When the argument of the Gamma function is a positive 
half-integer (given by the map z > n+ 1/2 > 0), the 
Gamma function simplifies to a double factorial: 


sa iad avila 


= [(2n — yna (5) jon 
= (2n — 1)!Vr/2". 


+ The product of all odd integers up to some odd positive integer 
n is often called the double factorial of n. It is denoted by n!!. 


Holder's Theorem 


A large number of definitions have been given for the 
Gamma function. Although they describe the same func- 
tion, it is not entirely straightforward to prove their equiv- 
alence. Instead of having to find a specialized proof for 
each formula, it would be highly desirable to have a gen- 
eral method of identifying the Gamma function given any 
particular form. 


One way to prove equivalence would be to find a differ- 
ential equation that characterizes the Gamma function. 
Most special functions in applied mathematics arise as 
solutions to differential equations, whose solutions are 
unique. However, the Gamma function does not appear 
to satisfy any simple differential equation. Otto Holder 
proved in 1887 that the Gamma function at least does 
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not satisfy any algebraic differential equation by show- 
ing that a solution to such an equation could not satisfy 
the Gamma function’s recurrence formula. This result is 
known as Holder’s theorem. 


A definite and generally applicable characterization of 
the Gamma function was not given until 1922. Harald 
Bohr and Johannes Mollerup then proved what is known 
as the Bohr-Mollerup theorem: that the Gamma function 
is the unique solution to the factorial recurrence relation 
that is positive and logarithmically convex for positive z 
and whose value at 1 is 1 (a function is logarithmically 
convex if its logarithm is convex). 


The Bohr-Mollerup theorem (discussed in the next chap- 
ter) is useful because it is relatively easy to prove loga- 
rithmic convexity for any of the different formulas used 
to define the Gamma function. Taking things further, 
instead of defining the Gamma function by any partic- 
ular formula, we can choose the conditions of the Bohr- 
Mollerup theorem as the definition, and then pick any 
formula we like that satisfies the conditions as a starting 
point for studying the Gamma function. 


The constant 7 = 3.14159... represents the ratio of the 
circumference of a circle to its diameter. The number 
m, like many other fundamental mathematical constants 
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such as e = 2.71828..., isa transcendental number. Both 
m and e possess an infinite number of digits which appear 
to have no orderly pattern to their arrangement. Tran- 
scendental numbers cannot be expressed as the root of 
any algebraic equation that contains only rational num- 
bers. In general, transcendental means nonalgebraic. 


To be clear, consider the numbers v2 and r. Irrational 
means ‘cannot be expressed as a fraction'. V2 is an ir- 
rational number but is not transcendental, whereas 7 is 
both transcendental and (is therefore) irrational. Any al- 
gebraic number raised to the power of an irrational num- 
ber is necessarily a transcendental number. 


A transcendental function does not satisfy any polyno- 
mial equation whose coefficients are themselves polyno- 
mials, in contrast to an algebraic function, which does 
satisfy such an equation. A transcendental function f(x) 
is not expressible as a finite combination of the algebraic 
operations (employing only rationals) of addition, sub- 
traction, multiplication, division, raising to a power, and 
extracting a root applied to x. Examples include the 
functions log, sing, cosg, e” and any functions con- 
taining them. Such functions are expressible in algebraic 
terms only as infinite series — in this sense a transcenden- 
tal function is a function that "transcends" algebra. 
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The following functions are transcendental: 


fila) = a” 

fot) Ss E701 
fsla) = =” 

Ly AS ae 

jolt): = loge, c#0,1 
jez) = sinz 


Formally, an analytic function f(z) of one real or complex 
variable z is transcendental if it is algebraically indepen- 
dent of that variable. This can be extended to functions 
of several variables. 


Hölder’s theorem states that the Gamma function does 
not satisfy any algebraic differential equation whose coef- 
ficients are rational functions. The result was first proved 
by Otto Hölder in 1887. The theorem also generalizes to 
the q-gamma function. 


Theorem 3.1. (Hélder’s Theorem). There exists no 
non-constant polynomial P(x; yo, Y1,---;Yn) such that 


P(z;T(z),T’(a),...,T™(x)) =0 


where Yo, Yı, -- - , Yn are functions of x, T(x) is the Gamma 
function, and P is a polynomial in yo, yı, ..., Yn with co- 
efficients drawn from the ring of polynomials in x. That 
is to say, 


P(g; Yo, Y1,--- Yn) E D OTE A (a0,a1,.-an) (©) f (yo) Panes (Yn) 
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where the (ao, a1,..., an) index all possible terms of the 
polynomial and the A(ag,ai,..,an)(©) are polynomials in x 
acting as coefficients of polynomial P. The Aao a...an) (£) 
may be constants or zero. 


Example 3.1. If P(x; yo, Y1, Y2) = 2°y2+ryi+(x?—a7) Yo 
then Ao,1) (2) = 2°, Avo,1,0) (x) = x and Aao ols) z 
(x? — a?) where a is a constant. All the other coefficients 
in the summation are zero. Then 


PRIT =a? f" + af' + (x? — 07) f =0 


is an algebraic differential equation which, in this exam- 
ple, has solutions f = J,(x) and f = Y,(zx), the Bessel 
functions of either the first or second kind. So 


P(x; Jo(x), J42), Jg(@)) = 0 


and therefore both J,(x) and Y,(«) are differentially al- 
gebraic (but are algebraically transcendental). 


The vast majority of the familiar special functions of 
mathematical physics are differentially algebraic. All al- 
gebraic combinations of differentially algebraic functions 
are also differentially algebraic. Also, all compositions 
of differentially algebraic functions are differentially al- 
gebraic. Holder’s theorem simply states that the gamma 
function, I(x), is not differentially algebraic and is, there- 
fore, a transcendentally transcendental (or equivalently, 
hypertranscendental) function. 
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A hypertranscendental function is a function which is not 
the solution of an algebraic differential equation with co- 
efficients in Z (the integers) and with algebraic initial 
conditions. The term was introduced by D. D. Morduhai- 
Boltovskoi in 1949. Hypertranscendental functions usu- 
ally arise as the solutions to functional equations, the 
Gamma function being one example. 


A functional equation is any equation that specifies a 
function in implicit form. Often, the equation relates the 
value of a function at some point with its values at other 
points. For instance, properties of functions can be de- 
termined by considering the types of functional equations 
they satisfy. The term functional equation usually refers 
to equations that cannot be simply reduced to algebraic 
equations. 


A few simple lemmas will greatly assist us in proving 
Holder’s theorem. 


In what follows, the minimal polynomial of an algebraic 
number ¢ is the unique irreducible monic polynomial of 
smallest degree m(x) with rational coefficients such that 
m(¢) = 0 and whose leading coefficient is 1. The min- 
imal polynomial of an algebraic number a divides any 
other polynomial with rational coefficients p(x) such that 
p(a) = 0. It follows that it has minimal degree among 
all polynomials f with this property. 
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The following table lists some algebraic numbers ¢ along- 
side their minimal polynomials m(x), computed using the 
Mathematica function MinimalPolynomial: 


Ç | m(x) 
2 g—2 
91/2 r2 
21/3 r3 
21/3 45 x3 — 15x? + 75x — 127 
21/2 4 31/2 | zf — 10x? +1 


Lemma 3.1. The minimal polynomial of an algebraic 
number Ç is unique. 


Proof. If we had two such polynomials, they must both 
have the same degree and the same leading coefficient 1, 
and so their difference is a polynomial of smaller degree 
which still gives 0 when applied to ¢. But this would 
contradict the minimality of m. 


Lemma 3.2. If p is some polynomial such that p(¢) = 0, 
then m divides p. 


Proof. By definition, deg(p) > deg(m). We may write 
p = qm-4r for some polynomials q,r, such that deg(r) < 
deg(m). Then since m(¢) = 0 and p(¢) = 0, we have 
that r(¢) = p(¢) — ¢(¢)m(¢) = 0, which contradicts the 
minimality of m, unless r(x) = 0. Therefore r(x) = 0 
and m divides p. 
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Proof. (Hélder’s Theorem). Assume the existence of 
P as described in the statement of the theorem, that is 


P(z;T(z),I’(x),...,T™(x)) =0 


with 
P(x; Yo, Yl,- Yn) a DET, Alaan (2) r (yo) D RS a (Yn). 


Also assume that P is of lowest possible degree. This 
means that all the coefficients A(ao,a1,...an) have no com- 
mon factor of the form (x — y) and so P is not divisible 
by any factor of (x — y). It also means that P is not the 
product of any two polynomials of lower degree. Consider 
the relations 


P(x +1;r(£ +1), r®(£+1),...,T®(£ + 1)) = 
= P (z + 1; 2T (2), [20 (2), Era), ..., [T(E ®) 
= P (x + 1; al (x£), sT® (x) +T(£),..., 2r” (x) + nre) (2)) 


so we can define a second polynomial Q, defined by the 
transformation 


Q(T; Yo Yrs -- Yn) = P (£ + 1; cy, TYL + Yo, +--+, TYn + NY(n-1) 


and Q (x;r(x), I" (£),...,T™®(x)) = 0 is also an alge- 
braic differential equation for T(x). Q and P both have 
the same degree and P must divide Q otherwise there 
would be a remainder and that would mean P was not of 
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minimal degree. Since (£yn + NYn—1)"" = gryte +., 
the highest degree term of Q is 


ghotiatethn Aa a hn)(@+1)- (yo) = (yr)... --(Yn)* 


where (ho, h1,...,n) are the exponents of the highest 
degree term of P. 


Let R(x) be the ratio between P and Q: 

Q(T; Yo. Yts «+ +5 Yn)= P (£ + 1; 2y0, cy1 + Yos- - - LYn + NYm-1)) 
= R(x)P(x; yo, Y1,- --, Yn)- 

Suppose R(x) has a zero, say y. Substitute y into 


P(y + 1; Yy0, YY + Yor +++ YYn + NYm-1)) = 0- P(T; Yo, Yi, +--+, Yn) = 0. 


This last equality indicates that (z—(y+1)) is a factor of 
P, contradicting the assumption that P was of minimal 
degree. 


Now consider the two leading terms, which must satisfy 
the equality 


R(x) Aho... ha) (£) (o ari (n) 
a GOO NG + 1) ` (yo)? REE (Yn) 
Ria) Ags. yn 8 A a E I: 


This equality cannot be satisfied for arbitrary x if R = c, 
where c is constant. Therefore, no such P exists and T (x) 
is not differentially algebraic. 


Bohr-Mullerup Theorem 


Theorem 1.3 shows that the Gamma function extends the 
factorial function from the set of natural numbers N to 
real non-null positive numbers. But the Gamma function 
is not the only way to do so. Consider, for instance, 
functions of the form cos(2max)I'(x), where m is any 
non-null integer. We already know T(x + 1) = zT (x), so 
it follows 


cos(2ma(a + 1))P(@ + 1) = xcos(2ma(x + 1))P (x) 
= x cos(2mrz + 2mr)I (x) 
= xcos(2mrx)T (x) 
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thereby satisfying the functional equation. Hadamard 
proposed another alternative: 


which gives F'(n+1) = n! for alln € N. In principle, there 
are an infinite number of possibilities, since we may draw 


any curve through the points (1, 1), (2, 2), (3, 6), (4, 24), ... 


and assume it represents some function that returns fac- 
torials at integer values. The purpose of this section is 
to explain why we consider the Gamma function as the 
function which extends the factorial function to all real 
x > 0. We seek a condition which implies that any such 
function cannot be anything else other than T(x). 


The functional equation does not only apply to the natu- 
ral numbers; the relation P(a+1) = zT (x) is valid for any 
x > 0. This is a stronger requirement than (n+ 1) =n! 
alone is, since it implies that the values of I in any range 
|x, x +1] determine the value of I on the entire real line. 
This is a fair amount of restriction, but not enough to 
meet our goal. 


To distinguish the Gamma function amongst all the pos- 
sible continuations of the factorial function, the notion of 
convexity is introduced. A standard definition of a con- 
vex function is a necessarily continuous function whose 
value at the midpoint of every interval in its domain does 
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not exceed the arithmetic mean of its values at the ends 
of the interval. Convex functions can be defined formally 
as follows: 


Definition 4.1. A function f : (a,b) > R is called con- 
vex if and only if for any x,y € (a,b), we have 


FAx + (1 Ajy) < Af (x) + (1 A) FY) 
for all A € (0,1). 


We have parametrized the interval (x,y) as {Ax + (1 — 
Ajy | 0 < à < 1}. The above definition means that 
as you move from x to y, the line joining (x, f(x)) to 
(y, f(y)) always lies above the graph of f. 


convex concave 


X 


Figure 4.1: Convex versus concave functions. 


Any convex function on an open interval must be contin- 
uous, which we will now show. 


Lemma 4.1. Any convex function f : (a,b) > R is 
continuous. 
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Proof. We shall follow the proof given by Rudin. The 
proof uses geometric language. See Figure 4.1 for the 
corresponding diagram. 


Supposea<s<x2<y<t<b. The points should not 
be considered fixed, just that the preceding inequalities 
are satisfied. Write S for the point (s, f(s)), and similarly 
for x, y and t. Convexity implies X must lie below SY, 
and Y is above the line SX. Also, Y is below the line 
XT. As y > xt, the point Y is sandwiched between 
these two lines, and hence f(y) > f(x). Left-hand limits 
are handled similarly, and continuity of f follows. 


As defined, convexity requires openess of the domain; 
otherwise, for instance, we may have a convex function 
on [a,b] which is not continuous, through the presence of 
a point discontinuity at a or b. 


Given two points x and y in the domain of a real function 
f, we may form the difference quotient n, which is 
equivalent to the slope of the line segment from (x, f()) 
to (y, f(y)). For convex functions, the difference quotient 
always increases as we increase x and y: 


Proposition 4.1. If f : (a,b) > R is conver and if 
a<s<t<u<b, then 


f(t) — f(s) < fu) — F(s) < Flu) = FE 


t-s u—s u—t 
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Proof. Leta <s <t< u< b. Then 
ettn f(s) + (1 Fw). 


Let A = “!, (When t = u, À = 0 and when t = s, A = 1.) 


Then 
astaw 2254 (1-22) 
u—s u—s 
_ (u=t)s+ (u-— s)u— (u — t)u 
7 u—s 


us — ts + u? — us — u? + ut 


u—s 


ut — ts 
= 


U— Ss 
It thus follows 


u—t 


FŒ) < fu) + 


(f(s) — f(u)). (4.1) 


U—S 


s—t 
u—s 


Use of the identity ut =1+ in equation 4.1 gives 


s—t 


FO) < f(s) + 


(f(s) — flu). (4.2) 


u—s 
Rearrangement of equation 4.1 gives, 


fu) — F) 


u—s u—t 


d 
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whereas rearrangement of equation 4.2 gives, 


f(t) — fs) — flu) — Fs) 


t-s a u—s 


Hence 


FŒ) — F(s) < f(u) — F(s) 


t—s u—s u—t 


as was meant to be shown. 


If we take the logarithm of both sides of the convex func- 
tion f(x) = x”, we get log f = 2logx. Now, log x is not a 
convex function (for each pair of points, the line segment 
joining them lies below the graph, not above it); so log x 
is in fact concave, not convex. Thus, log f is not convex 
either. But if we consider h(x) = e”, which is convex, we 
have that logh = x, which is also convex. We use the 
terminology log-convex to describe h, which is a stronger 
property than convexity is. To show this, we need the 
following preliminary lemma: 


Lemma 4.2. Any increasing convex function of a convex 
function is convex. 


Proof. Let f : (a,b) — (h,k) be convex and let g : 
(h,k) — R be convex and increasing, i.e., £ < y > 
g(x) < gly). By the convexity of f, for anya < z < y < b 
we have 


f(Ar + (1 — Ay) < Af(£)+ (1 - Af (y) 
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for any A € [0,1]. Furthermore by the convexity of g we 
get, 
g(f (At + (1 = A)y)) < 9AF (2) + 1 - A) FY) 
< Ag(f(x)) + (1 — A)g(F(y)). 


Hence g o f is convex. 


T(x) Log(I'(x)) 


6 6 
5 5 
4 4 
3 3 
2 2 
1 1 


Figure 4.2: The graphs of T(x) and logI'(x) plotted over (0,4). 


Proposition 4.2. Given a function f : (a,b) > R, if 
log f is convex, then so is f itself. 


Proof. If log f is convex, then by lemma 4.2, since e” is 
increasing and convex, we have that elf = f is also 
convex. 


We now return to the Gamma function. It is obvious 
from figure 4.2 that T(x) with x € (0,00) is convex. In 
fact, T(x) increases steeply enough as x > oo that r : 
(0,00) => R is log-convex, which we shall momentarily 
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prove. To prove that logT is convex (which implies that 
T is convex) we need an inequality known as Holder’s 
inequality, which we prove in kind. To this end, consider 
the following: 


Lemma 4.3. Jfa >0,b<1,5>0,t>0ands+t=1, 
then 
(1-a) (1 — b) < 1-— at. (4.3) 


Proof. Consider f(a) = (1 — a) (1 — 6)’ + ab where b is 
fixed. Then 


f'(a) = —s(1 — a)®71(1 — b)! + sa? 28! 


from which it follows 


1 1-s 1 t 
=O => (<-1) = (5-1) > a=b. 
Now 
f"(a) = sls — 1) =a)? "= b) + s(s = 1a? 7b 
from which we see 
Fal =s(s-1) [1 -— a)! + | <0. 
This means f attains its maximum at a = b, and we have 


f(a)|,__,=l-a+a=1. 


Hence f(a) < 1 which gives inequality 4.3. 
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Lemma 4.4. (Holder’s inequality). Let p and q be 


positive real numbers such that ++ 7 = 1. Then for any 


integrable functions f,g : [a,b] + R, we have 


< (f ira) (firar) 


| [ tosa 


1/4 
(4.4) 


Proof. s + ; = 1 > (1 — p)(1 — q) = 1, thus inequality 
4.3 is satisfied. Now, using inequality 4.3 as well as the 
Cauchy-Schwarz inequality, we have in the limit € > OT, 


b— 


< f:g dz 


[tou 


+ | f()9(b) |e 


a 


1/p 1/q 
< ( gipas) (12wirar) +|f(0)| -|g(d)| -€ 
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-(f fP dx |f) J cc glt dx — [g(0)|"- ) ' 


+ |f(d)| - |g(b)| -€ 

_ \/p in 
(2 ra) (stores) (1-a) (i-t) 
+ | f(d)| s®l € 


1/q 
b |fllglļe 
< ( S, \f|? 1) He tar] (: mise] 


and the proof is complete. 


+ IFO) - |g(b)| -€ 


-( [ rar) “( / igrae) k 


Now consider that 


Theorem 4.3. T : (0,00) > R is log-convez. 


Proof. Let 1 < p < œ and EE = 1. As in Rudin, 


41 


consider 


(ste) 


[x «/p+y/q—-1) o~t qt 
0 


E pt/Ppy/a¢—1/Pp-1/4e—t/Pe—t/4 dt 
0 


(= 1 e~) 1/p let 1/4 qt 
l wre) 


es 1/p 3% 1/q 
<( i Plet a) ( | tolet 7 
0 0 


= P(e)" Py)". 


IA 


Let À = a and hence 1 — à = = Then à € (0,1), and 


(Ag + (1-A)y) < P(e Ty)? 
logT (Az + (1 — A)y) < log [F(£)^ Py) ] 


for any x,y € (0,00). Hence logIT is convex. O 


A famous theorem of Bohr and Mullerup says that propo- 
sition 1.2, theorem 1.3 and theorem 4.3, taken together, 
distinguish T as the function which extends the factorials 
to all real x > 0. 


Theorem 4.4. (Bohr-Mullerup). If f : (0,00) > 
(0,00) satisfies 
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2. f(x +1) = zf(x), and 
3. log f is convex, 
then f(x) =T(x) for all x € (0,00). 


Proof. Rudin’s proof is very elegant. Since we have al- 
ready shown T to satisfy conditions 1 through 3, it suf- 
fices to prove that f(x) is uniquely determined by these 
conditions. Furthermore, condition 2 asserts it is enough 
to prove this only for x € (0,1). 


Set y = log f. Condition 1 says that y(1) = 0 and 
condition 2 says 


plz +1) = v(x) + loge. (4.5) 


Condition 3 means that y is convex. 


Let 0 < x < 1, and let n € N. Consider the difference 
quotients of y; let s=n,t=n+1landu=n+1+2 to 
get 

p(n+1+2)—y(n+1) 


y(n+1)—Yy(n) < - l 


The difference quotients of y with s =n+1,t=n+1+2 
and u = n + 2 give 


yp(nt+1+2)—-y(n+1) 


< y(n + 2) — p(n +1). 


Note that by equation 4.5, p(n+1)—y(n) = logn, and in 
addition by condition 2, y(n +2)— y(n +1) = log(n+1). 
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Combining these gives 


yp(nt+1+2)—-y(n+1) 


logn < 


< log(n + 1). 


Repeatedly applying equation 4.5 gives 


p(z+n-+1) 
= p(x +n) + log(z +n) 
= y(x +n — 1) + log(z +n) + log(x +n — 1) 
= p(x +n — 1) + log|( 
(cx +n— 2) + logl(a +n)(z +n -— 1)(z +n -— 2)] 


= p(x) + log|(x + n)(x +n — 1): +- (x + 1)a]. 
Also by equation 4.5, we have y(n + 1) = log(n!). So 


yp(nt+1+2)—-—y(n+1) 


~ lela) + log|(x + n) --- (z + 1)a] — log(n!)] 


giving us 
logn < 4[p(x) + log[(z +n): -- (a + 1)a] — log(n!)] < log(n + 1). 
Multiplying through by x yields 


logn” < y(x) + log|(x +n)--- (x + 1)a] — log(n!) < log(n + 1)’. 
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Subtracting log n” from each term gives 


0 < y(x) + log[(x + n)--- (a + 1)2] — log(n!) — logn” < log(n + 1)” — log n”. 


1 
< zlog (+5). 
n 


Now let n — oo, so that log (1 + +) — 0, and hence 


Simplifying, 
n'n” 


0 < y(x) — log z(z +1) (£ +n) 


n'n? 


a E 


n—> o0 


In any case y is uniquely determined and the proof is 
complete. 


The last equation in the above proof brings us to an al- 
ternative definition for the Gamma function. In another 
(earlier) letter written October 13, 1729 also to his friend 
Goldbach, Euler gave the following equivalent definition 
for the Gamma function: 


Definition 4.2. Let 0 < x < o and define 


; n'n? : n'n? 
D(a) = Mni = = limn +00 [| : 
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This approach, using an infinite product, was also chosen 
by Gauss, in 1811, in his study of the Gamma function. 
Using this formulation is often more convenient in estab- 
lishing new properties of the Gamma function. 


Proof. Since we are dealing with a definition, this proof 
is only to show consistency with the previous definitions 
of I already given. The proof of theorem 4.4 shows that 
n'n? 

g(x +1)--- (z +n) 


p(x) = lim log 


noo 


for 0 < x < 1. Since the log function is continuous, we 
can exchange lim and log this way: 


n'n? 


u(x +1) (x+n) 


g(x) = log lim 


N—>Co 


Exponentiating both sides gives 


nin” 


x(a +1)---(a+n) 


T(x)= lim 


n—> o0 


(4.6) 


for 0 < x < 1. Equation 4.6 also holds for x = 1, in 
which case we have 
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Thus, equation 4.6 holds for 0 < x < 1. Using proposi- 
tion 1.2, we see that 


nin 
T 1)=x<x li 
Ca | GaGa) 
. z+n+1 nine 
= li 
nao n (z +1) (z +n) +n+1) 
Ear 14a]. jj nine 
= im 1+ 7 im e 
ninttt 
= lim f 
n> | (x +1) (x+n) +n+4+1) 


From this we see that equation 4.6 holds for 1 < x < 2 
as well; repeatedly applying this procedure shows that it 
applies for all x > 0, as required. 


The Beta Function 


The integral in the definition of the Gamma function 
(definition 1.1) is known as Euler’s second integral. Now, 
Euler’s first integral (1730) is another integral related to 
the Gamma function, which he also proposed: 


Definition 5.1. Beta Function. For Re(x), Re(y) > 0, 
define 


B(a, y) := age — t)” dt 
=2 u sin(t)?*~' cos(t)®t dt 
P(x) (y) 


Fa 
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This integral is commonly known as the Beta function. 
The definition above involves three equivalent identities — 
an integral over trigonometric functions, an integral over 
polynomials and a ratio of Gamma functions. The Beta 
function is symmetric. 


To interrelate these various representations of the Beta 
function, consider the product of two factorials written 
in terms of the Gamma function: 


m!n! -|/ eu” au | e "uv" dv. 
0 0 
Let u = x? and v = y? so that 
m!n! = af en gmt ar | eV yt] dy 
0 
a fe (x2?+y?) io gh dz dy. 


Switch to polar coordinates with x = rcos@ and y = 
rsinð, 


OO 
a2 i 
m!n! = f I e" |r cos 6/?"t" |r sin 0|” r dr dO 
o Jo 
oo 5 2T 
= il erT as ar f | cos?” t! @ sin?"t! 0| d0 
0 0 


oo T/2 
= af en pAmtnt)+1 ar f cos2™+! @ sin2"+! 6 dé. 
0 0 
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Now make the substitutions t = r° and dt = 2r dr in the 
first integral, 


oo T/2 
min! =2 f eO iE ar f cos?”+! @ sin?” tt 6 da 
0 0 
a /2 
=2-T(m+n-+2) f cos?™ +i 9 sin?” t! 9 d0 
0 
a /2 
=2(m+n +1)! i cos" 8 sine? Ode. 
0 
In terms of the Beta function we then have 


T/2 
B(m+1,n+1)= 2 | cos?™ +! @ sin?" 6 dO 
0 


m! n! 
~~ (m+n++1)! 


Adjusting the arguments then gives 


From the trigonometric form 


a /2 
Becig@ => i sin(t)2"+? cos(t)2#+! dt 
0 
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we recover the integral over polynomials by making the 
change of variables u = cos? t and du = 2costsint dt. 


Doing this, 


Be +1.y+1)= | (1 — u)” (u)” du 


so readjustment of the arguments and relabeling gives us 


Boae | (1 — 07-1) dt. 


Directly from the definition follows the beta function func- 
tional equation 


z 
B(x +1,y) = pa 


Proof. 


— P&r) — arau) _ _« 
Bia +1,9) = Tent) = erret) 7 sty” (29). 


By substituting z = y = 5 into the definition, we obtain 
the following: 


Corollary. ['($) = yr. 
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Proof. Considering 


P(x)I (y) 


T/2 
DEE) sin 6)?! (cos 6)74~! dé 
Te I (sin 8)2*-(cos 8) 


the special case 7 = y = z gives 


CO 


Since I is positive for all x € (0, 00), we have ['(3) = y7 
as required. 


Using this result, we can easily perform an integral very 
fundamental in probability theory: 


Proposition 5.1. 


f e” dr = Vr. 


oO 


Proof. We have T(4) = yr. By definition, 


1 CO 
T G) = f te~ dt. 
2 0 


Substituting t = x?, dt = 2x dx yields 


1 i 2 
r(5) = f g te ™ Qaeda = yT. 
0 
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Therefore m 
T 
f e`” dr = vr 
0 2 
Since e7?" is symmetric about x = 0, we have 


| is] e dr = V7. 
2 0 


OO 


Wallis’s Integrals 


The following integrals (Wallis’s integrals) 
T2 T/2 
W,, = f sin” 8 d0 = f cos” 0 d8 
0 0 


may be computed by means of the Beta and Gamma 
functions. By the definition of the Beta function, we 


have 
1 n+1 1 
Wa = 53B —=— >53 
2 ( 2 5) 
which gives rise to the two cases n = 2p + 1 and n = 2p. 
For the odd values of the argument n: 


r(p+1)r(1/2) sp! (1/2) 
2 


1 
Wəp+1 = 5p(» rl, 2T (p + 3/2) (2p + 1)F(p + 1/2) 


53 
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Using the formula 


T(n + 1/2) = 


produces the result 


W. 2P p! 4Pp!? 
Pete ae es ar pari)! 


For the even values of the argument n: 


1 1 1\T(p +1/2)0 (1/2) 
Way = 3B(p3 5) T(p +1) 
and 
Woes 1-3-5-(Qp-1) — (Qp)!a 


2P+1p! Appl? 2° 


We obviously have 


Wr42 = 2B (So 3) E ew, = Gag) Wn 


according to the Beta function functional equation. 


Note that 


1 a+l 1 
Wa = =B = 
2 ( 2 3) 
works for any real œ > —1 and we can therefore deduce 
using the definition of the Beta function (respectively 
with a = —1/2 and a = 1/2) that 
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le 1 2dt — T(1/4) 
sin o Vl—-# 2/27 


T/2 1 2% dt (2r)? 
v sin 0 d = = : 
0 oa bee  T?(1/4) 


The product of those two integrals allows us to derive a 
relation due to Euler: 


: l Paka i 
o vVl=tiJo VI-# 4 


Wallis’s Product 


Let’s establish an infinite product for 7/2 known as “Wal- 
lis’s product.” 


Theorem 7.1. 


nm 224466 2k 22k 
2 133557 2k—12k+1 


By this is meant that if P, is the product of the first n 
factors on the right-hand side, 


lim P, = a3 
2 
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Proof. Directly from the definition of the Beta function, 


11 ee r 
B(n+5.5) : J sin?” z dz = VALUR UA 
0 


2°92 2(nl) 
1 me Jan! 
B 1,=): int gde = — 
(n+ 5) is sin x dx 2T(n + 3/2) 
where n = 0,1,.... (7.1) 


Hence, the quotient of these two integrals is 


Tes sin gde — T(n+3)T(n +8) 
fe sint edz n! n! 
— 2mm+1l2n-12n-1 3317 


Im WwW M-22 4222 
Ir 


TBR? 


(7.2) 


We shall now show that the left-hand side of equation 7.2 
approaches 1 as n — oo. By equation 7.1 formed for n 
and for n — 1 we have 


Since 0 < sinx < 1 in the interval (0,7/2), we have 


Ja 2. 2, ogee 
0< Wa sin?” t! zdz < ne sin?” zdz < z/ sin?! gdz. 
0 0 0 
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Dividing this inequality by the first of its integrals and 
allowing n to become infinite, we have by equation 7.3 
that the left-hand side of equation 7.2 approaches 1. 


Hence, 
n= 2 
Also 7 j 
+ T 

lim Pons = lim 2 Pon = 

n—0o n=>œ 2n +1 2 
and the proof is complete. 
Corollary. 

(n!)?2?% 


ii (2n)! y/n SAN 


Proof. To prove this, multiply and divide the right-hand 


side of the equation 


2 22 Qn Nn 
“13. OA TH Se 


Pop, 


by 2:2---2n-2n, thus introducing factorials in the denom- 
inator. If then factors 2 are segregated in the numerator, 
the result becomes apparent. 


Product & Reflection 
ormulas 


Using the fact that (1—t/n)” converges to e™* as n — on, 
one may write 
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Integrating by parts, we have 
1 n 
T(z) = lim —.- zf t(n — t dt 
0 


1 ies n 
Sfin D f ptn- dt 
no n” z(z + 1)--- (z +n- 1) 0 


o n 1 2 n 
= lim ee d 
tim = (5) (5) (=) 


To evaluate the above limit, we insert convergence factors 


e~7/* to get 
1 zim ree TT +5 e77/k 
T(z) n—0o 
k=1 
= lim er (lt1/2+--+1/n—log n) eo 2/k 
nN—- Oo IT( ; 


We shall shortly prove that 1 + z +-+ 4 — logn ap- 
proaches a positive limit y, known as the Euler constant, 


so that 3 
=a ze” I] 142 e (8.1) 
T(z) k i l 


k=1 
which is Weierstrass’ product form of the Gamma 
function. 
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Euler’s constant, also known as the Euler-Mascheroni 
constant, has the numerical value 


y ~ 0.577 215 664 901 532 860 606512.... 


Using the Weierstrass identity to define an extension of 
the Gamma function to the left half-plane, we get 


1 oT 2 sin TZ 
a neg, eee 
uS a I ( =| ap 


where we used the identity 


which is proved in Appendix A. 


Thus 


T(z)P(—z) = (8.2) 


Zsin Tz 
and given that [(1 — z) = —zI'(—z), 


irra bea 


(8.3) 


sin TZ 


which is Euler’s reflection formula. 


Two immediate results which follow from equation 8.3 
are 
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1. T is zero-free, 


2. [(1/2) = yr. Applying T(z + 1) = zI(z), we also 
find T(3/2) = 4yr, r(5/2) = 37/4, and so on. 


Lemma 8.1. (Euler’s constant y). Ifs, = 1+ 4+ 
-+ 1 — logn, then liMnp—oo Sn exists. 


Proof. tn = 1+4 +--+- +1/(n— 1) — logn increases with 
n since, in the geometric sense, t, represents the area of 
n — 1 regions between an upper Reimann sum and the 
exact value of {/"(1/x) dz. We may write 


k=1 
and 
=} 1 
lim tn = ——] 1+-— 
k=1 
Now since 
r? r q’ 
log(1 == = Pree 
og(l+z)=z2 ae ee 
the series above converges to a positive constant since 
1 1 1 1 1 1 
`] {a = eS | Saye es, 
is ae ( :) Dk? 33" Akt = OR? 


This proves the lemma, because limp so Sn = liMp—so0 tn- 
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Using equation 8.3, the gamma function I(r) of a rational 
number r can be reduced. For instance, 


r(5) = ATG 


Half-Integer Values 


From the functional equation [(n+1) = nI (n) the value 
of the Gamma function at half-integer values is deter- 
mined by a single one of them; one has 


from which it follows, through repeated application of 
the functional equation, that for n € N, 


r(5 +n) -v -aD y 


2N 
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Similarly, 


n\ (n-2)! 
r(3) = CDI yr. 


10 


Digamma and Polygamma 
Functions 


An entire theory revolves around Digamma and Polygamma 


functions, which we shall not pursue in great depth here. 
We briefly define them, then first show their involvement 
in the derivation of an expression for the derivative of the 
Gamma, function which follows from Weierstrass’ prod- 
uct. Series expansions shall prove useful. 


The Digamma function is given by the logarithmic deriva- 
tive of the Gamma function: 


d T(z) 


p(z) = z, os Fe) = 


The nth derivative of Y (z) is called the polygamma func- 
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tion, denoted w,(z). The notation 


polz) = V(z) 


is often used for the digamma function itself. 


Begin with the Weierstrass form 


Lig = 2e II (: + ) i“ ; 


Take the logarithm of both sides, 
E S log tE Sn 
k k 


k=1 


Now differentiate, 
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It immediately follows that 


-rw (141-55) 


Given the expression for I”(z), the minimum value £m of 
T(x) for real positive x = £m occurs when 

I'(£m) = P(@m)bo(@m) = 0. 
In other words, when 


Numerically, this is solved to give £m = 1.46163.... 
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We have that 


oibie 


k=1 ) 
=y z #0, -—1, 
n a) 


(10.1) 


-7 4 = (5 Hehe S z#0,-1,-2,... 


If we differentiate relation 10.1 many times, we find 


v(e) = PI") -Tr e) _ s5 1 


T(z) = k+z—-1) 
(10.2) 
Uw" (z =-)'_——_,, 
= (k+2—1)8 
in, ve (Pha! 
miels CE (10.3) 
k=1 


where the Y, = Y™ functions are the polygamma func- 
tions, 


v(e) = 2 (loa(t(2)) 
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The series expansion relation 10.1 suggests 


vern-we)= 3° (mh) 


k=1 


which gives the recurrence formula 


1 
W(iz4+1) = V(z)+- 

z 

1 1 1 
R4 = Ü l per >1 

(A (2) z z+l1 z+n-1 Ga 
and by differentiating the first of these relations, 
B (—1)"n! 


Now, the Riemann zeta function, ¢(s), is a function of a 
complex variable s that analytically continues the sum of 
the infinite series 


The Riemann zeta function is defined as the analytic con- 
tinuation of the function defined for Re(s) > 1 by the sum 
of the given series. 


The values of the zeta function obtained from integral 
arguments are called zeta constants. The following are 
some common values of the Riemann zeta function: 
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1 
jas 
o=- 
1 1 
Hariss — 
¢(1) art OO 
1 1 r? 
2—1 = 
1 1 
68) =14 5+ yt 1.202... 
1 1 ni 
4)=1 nru x 1.0823... 
¢(4) 24 34 90 


Let’s consider the Basel problem in greater detail. The 
Basel problem is a famous problem in number theory, 
first posed by Pietro Mengoli in 1644, and solved by Eu- 
ler in 1735. Seeing that the problem had withstood the 
attacks of the leading mathematicians of the day, Euler’s 
accomplishment brought him immediate fame when he 
was twenty-eight. Over time, Euler generalized the prob- 
lem considerably, and his ideas were taken up later by 
Riemann when he defined his zeta function and proved 
its fundamental properties. The problem is named after 
Basel, hometown of Euler as well as of the Bernoulli fam- 
ily, who unsuccessfully attacked the problem. 


The Basel problem asks for the precise summation of the 
reciprocals of the squares of the natural numbers, i.e., a 
closed form solution to (with proof) of the infinite series: 
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aot i Ka D, d 
e — im t fee i ` 
A n2 n>+00 \ 12 22 n2 


Euler’s arguments were based on manipulations that were 
not justified at the time, and it was not until 1741 that 
he was able to produce a truly rigorous proof. 

Euler’s “derivation” of the value Li is clever and origi- 
nal. He extended observations about finite polynomials 
and assumed that these properties hold true for infinite 
series. By simply obtaining the correct value, he was able 
to verify it numerically against partial sums of the series. 
The agreement he observed gave him sufficient confidence 
to announce his result to the mathematical community. 


To follow Euler’s argument, consider the Taylor series 
expansion of the sine function 


sinft)_,_2 tae 
x 3! ! 


Now, the zeros of sin(x)/x occur at x = n -m where 
n = +1, +2, +3,.... Assume we can express this infinite 
series as a product of linear factors given by its zeros, 
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just as is commonly done for finite polynomials: 


sin(x) 


"S tre) 


IAA 


If we formally multiply out this product and collect all 
the x? terms, we see that the x? coefficient of sin (x)/x is 


1 1 1 C&A i 
= | | Desi In = 
(2 l 4r2 l 972 l ) =D 


n=1 


But from the original infinite series expansion of sin()/z, 
the coefficient of x? is —1/(3!) = —1/6. These two coef- 
ficients must be equal. Thus, 


1 ig 
ea oO 


n=l, 


Multiplying through both sides of this equation by —7? 
gives the sum of the reciprocals of the positive square 
integers, 

Zail, ET? 

ne Le 

Now, returning to the Polygamma function, from the re- 
lations 10.1 and 10.3 we have 


TT 


WSE 
Y1 (1) = ¢(2) = 7/6 
Wa(1) = —2¢(3) 
W,(1) = (-1)"*?n! C(n + 1) (10.5) 


Using the recurrence relation 10.4 allows us to compute 
these values for any positive integer, 


V(n) = i Ser D (10.6) 
= —y + Hn-1- 


The following series expansions are obvious consequences 
of relations 10.5 and of the series 


1 lo) 
eee — 1) k-1 
1+ > a 


Theorem 10.1. (Digamma series). 


T(1 +g) = -y + X (-1)*¢(k)z** |z| <1, 
= (10.7) 
(1+2) = ———-(y-1) + (1) (¢(k)—1)z"! |z| <1. 
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Series Expansions 


Finding series expansions for the Gamma function is now 
a direct consequence of the series expansions for the Digamma 
function: 


Theorem 11.1. For |z| < 1, 


log(T (1 + x)) = —yz 4 F ZA 


(11.1) 
log(I'(1 + x))= —log(1 + 2) — (y— La + ZZ, CECO- pt 
(11.2) 


Proof. Use term by term integration of the Taylor series 
10.7 and 10.8. 
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It follows easily from equation 11.1 and the functional 
equation that 


SER Se 
Tz) = zexp -5 ie ) (11.3) 


k=2 


Series formulas involving Ç(k) can be derived from equa- 
tion 11.1. For example, setting x = 1 gives 


log(P(2)) = -7+ a 


thus we arrive at a result due to Euler, 


n (DCC) 
k= k ; 


2 


y= 
Setting x = 1/2 into equation 11.1 yields 


i (r(3)) E crew) a 


so we have 
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Euler-Mascheroni Integrals 


Using the integral representation of I’(x) gives the inter- 
esting integral formula for Euler’s constant 


and from 


yy Tare -r & 1 
Eis T2(z) ny 


comes the relation 


P(r?) +r?) = r0") 


hence 
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We may continue this, computing the Euler-Mascheroni 
integrals — 


1 
T®(1) = -7 — Sn?y - 2¢(8) 
3 
PO (1) =o" + 04? + 86(8)7 + aot 


PON) = =a? = $4 = 206 (8)? — Fay = 246(5) — FC) 
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Duplication & Multiplication 
Formulas 


Theorem 13.1. (Gauss Multiplication Formula). 


o (24 22) = (en) DPT (2). 
y ye (2+) c(t =) (27) 0D pV a 


A brilliant proof of the multiplication formula was pro- 
duced by Liouville in 1855. We here present a modernized 
version of that proof. 


Proof. The product of Gamma functions on the left-hand 
side of equation 13.1 can be written 


S eet dt, fies pe aa dtz- -- i erage byn-1 dtp. 
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Rearranging, 


JE ibe at di en (tattatetin) gel? +(1/n)-1 pet ((n—1)/n)-1 dorat. 


Next, introduce a change of variables: 


The Jacobian is thus 


n—-1 


nq 
totz- tn 


so the integral can be written 


Eea] (1 Hts cg a Te =] 


z—1 
x( q” ) gtn) E A 1)/n)—1 ng” ~ dqdtz-- - dtp. 


ta--ty tot3:--tn 


Set s=totitg+---t+tr+q"/(tots---tn), so we have 


n ll Vee ao (emt EAS (1/n) aa Als gl 1)/n) 1 dqdty ++ dtn. 
(13.2) 
Now evaluate 


D af ee dtz dt3--- dt, 
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Obviously, 
d7 = k/n)—1 dt2 +- dtn 
ae, ee Te ee 
Now introduce a second change of variables, 


0" Mi tn), t3 = t3, ... tn = tn 


and 


S1 = t3 } t4 pesee] tn ti } q / Cta- trti). 


Now the Jacobian is 


E tts nae 
and we have 


d7 CO CO es 
da: = nq?! IA A Jo (at = 


1 
Senf aa 2 A Grae A 


k+1 


(1/n)— n-1 
plhin) —1 dt, dtz: dtn 
k+1 q” /t1 


k=2 


= —nZ. 
Thus, 
T= CE", 


To find C, we set q = 0 in the Z integral and in the 
equation above, then set them equal to each other to get 


(Eafe 
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Next we factor the integer n in the following way: 


n—1 
n= Jfa — expl2ari/n)) 
a=1 
2 — 1 
ED i ti i e. (13.3) 
n n n 


Here we used the identities 
2 2 
1 — exp (22m) a= cos an ¿sin = 
n 


sin 2i sin cos 
n n n 


n 
. r| . an, ar 
= 2 sin — | sin — — i cos — 
n n n 
and we also have 


- aT in Xo 
I sin" —icos—~| = (-i)* 'exp| —S a 
_ n n 


so that equation 13.3 follows. 


Now by the reflection formula 
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so by letting b in r(b) run from 1 to n — 1, we essen- 
tially perform the product given by C twice, so that 
CaO en and T = On) Og ie, 


Substitution in equation 13.2 gives 
T(z)P(z + 1/n)---T(zt+ (n — 1)/n) 
as nanyene f a gta dq 
0 


= n!/2=72 (27) "FT (n2) 


which completes the proof. 


Legendre’s duplication formula is found simply by setting 
n = 2 in the multiplication formula. 
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The Gamma and Zeta 
Function Relationship 


We have already shown in a previous section that 


1 ~ (=1)¥¢(k) 24 
res mas (x bp i ) 


k=2 


Now we will derive two more relationships between the 
Gamma function and the zeta function. 


Theorem 14.1. (Gamma and zeta functions rela- 
tion). 


oo 471 


C(z)T(z) = j aad dt for Re(z)>1. (14.1) 
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Proof. The integral definition of the Gamma function 


r= f tte" dt, 
0 


together with the change of variables t = ku (with k a 
positive integer) gives 


T(z) = (hu) eM du =H | eer du. 
0 0 


Rewrite this in the form 


oo 1 1 ik X oo = 
J — = u? > e~) du 
k=1 ke T(z) Jo 2l 
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For z = 2, this becomes 
2 co t 
=- dt. 
6 0 et — 1 


There is another important functional equation between 
the two functions, the Riemann zeta function functional 
equation: 


Theorem 14.2. (Functional Equation). The func- 
tion ¢(s) is regular for all values of s except s = 1, where 
there is a simple pole with residue 1. It satisfies the func- 
tional equation 


C(s) = 280! sin (S)ra —s)¢(l—s). (14.2) 


Proof. This can be proved in many different ways. Here 
we give a proof depending on the following summation 
formula. 


Let (x) be any function with a continuous derivative 
in the interval [a,b]. Then, if [x] denotes the greatest 
integer not exceeding x, 
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Since the formula is additive with respect to the interval 
(a, b], we may suppose n < a < b < n+1 so that we have 


fie —n—1/2)d'(x) dx = (b — n — 1/2) (0) — (a — n — 1/2)¢(a) — Ta d(x) da. 


The right hand side of equation 14.3 thus reduces to 
([b] — n)d(b). This vanishes unless b = n + 1, in which 
case it is d(n + 1). 


Let (n) = n~*, where s 4 1, and let a and b be positive 
integers. Then 


(14.4) 
Take the half-plane o > 1, a = 1, and make b > on. 
Adding 1 to each side, we obtain 


¢(s) = f 2 — 2 dz + : | = (14.5) 


s—-l 2 


Since |x] — x+ 1/2 is bounded, this integral is convergent 
for o > 0, and uniformly convergent in any finite region 
to the right of ø = 0. It therefore defines an analytic 
function of s, regular for o > 0. The right-hand side 
therefore provides the analytic continuation of ¢(s) up to 
o = 0, and there is clearly a simple pole at s = 1 with 
residue 1. 
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For 0 < g < 1 we have 


S [ dz 1 
2), xst 9 
and equation 14.5 may be written 
¢(s) = >| Ble da (Oe =) (14.6) 
0 


stl 


Actually equation 14.5 gives the analytic continuation of 
¢(s) for o > 1, for if we have 


rskl- Aes f Fod 


then fı(x) is also bounded, since, as is easily seen, 


k+1 


f(y) dy = 0 
k 


for any integer k. Hence 


i i 


gstl 


gst + (s +1) gst2 dx 


T1 Tı 


ae e "2 f(x) 


which tends to 0 as 71 > œ, £2 > œ, if o > —1. Hence 
the integral in equation 14.5 is convergent for o > —1. 
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Also it is easily verified that 


‘{e]-a2 +34 1 1 
= l dys > < 0). 
sf e A 


Hence 


C(s) =0f Hota (-l<o<0). (147) 


gstl 


Now we have the Fourier series 


1 sin NTT 
as pee ` ae 14.8 
[x] i 2 NT ( ) 


n=l 


where x is not an integer. Substituting in equation 14.7, 
and integrating term by term, we obtain 


sin2nne 4 
= 25 PS 
s (2n7)* sin y 
= D n [ ysti dy 


n=1 


4 


ST 


(2ny"{T(-s)}sin (FCs 


A | % 


i.e., equation 14.2. This is valid primarily for —1 < ø < 
0. Here, however, the right-hand side is analytic for all 
values of s such that ø < 0. It therefore provides the 
analytic continuation of ¢(s) over the remainder of the 
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plane, and there are no singularities other than the pole 
already encountered at s = 1. 


We have still to justify the term-by-term integration. 
Since the series 14.8 is boundedly convergent, term-by- 
term integration over any finite range is permissible. It 
is therefore sufficient to prove that 


sin NTE 4 
im J5 f o =O aN 


A—> o0 


Now 


a SE a _ [- cos dn] ™ s+1 fœ cos2nmx q 
A = 2nr JA est 


etl 2nngzs+i ~~ 


= O (sy) +O (gS z) = O (em) 


and the desired result clearly follows. 


Corollary. (Functional Equation). Let 


Ag) =m r(3) ets) 
be an analytic function except at poles 0 and 1, then 
A(s) = A(1— s). 
Proof. Changing s into 1 — s, the functional equation is 


ST 


(t= s) = 21r cos (F)ree(s). (149) 
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It may also be written 


G(s) = x(s)¢(1 — s) (14.10) 
where 


sl (1/2 — 8/2) 


x(s) = 2°! sin (F)ra _s)=n A 


and 
x(s)x(1-s)=1. 


The corollary is at once verified from 14.9 and 14.10. 
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Stirling’s Formula 


In this section we study how the Gamma function be- 
haves when the argument x becomes large. If we restrict 
the argument x to integral values n, the following result, 
due to James Stirling and Abraham de Moivre is quite 
famous and extremely important: 


Theorem 15.1. (Stirling’s formula). If the integer 
n tends to infinity, we have the asymptotic formula 


Tin+1) =n! ~ V2ann"e™. (15.1) 


For greater clarity in the proof of equation 15.1, we in- 
troduce several simple lemmas. 


Lemma 15.1. 


i pei a 1,2 
— = 10 T n= pyres 
z £ $ n 2n +1 
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Proof. Since the curve y = 1/z is convex, the area under 
the curve from x = n to x = n+ 1 is greater than the 
area of the trapezoid bounded by these two ordinates, the 
x-axis, and the lower tangent to the curve at the point 
(n+ 1/2,2/(2n+1)), i.e., the lower Riemann sum. Note 
that the area of a trapezoid is equal the product of the 
length of the base by the length of the median, which, in 
this case, is 


I a: -2 
n+ Ons 
Lemma 15.2. 
n! 
lim a, = lim ———— erists. 


Proof. We have 


An+1 e 


since by lemma 15.1, 


1 1 
(n+5) tog (1+) >1. 
2 n 


Since a, > 0 for all n, the proof is complete. 
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Lemma 15.3. 

lim a, > 0. 

n—> o0 
Proof. Use an argument similar to that of lemma 15.2, 
but this time use an upper Riemann sum and compare 
areas with the function y = log x. We break the interval 
up by bounding a series of trapezoids by two rectangles, 
one at each end of the interval. The altitudes of the two 
rectangles at the endpoints are 2 and logn. Thus, the 
area of the trapezoids and the two rectangles is 


1 
1 +log2+log3+---4 log(n — 1) + 5 logn 
= 1 + logn! — log yn 


The area under the curve is 


f logx dz = nlogn = n + 1 = log(n/e)” + 1. 
1 


Hence, 
log | Kai i 
Ee oa 
ee ee ee 
n! 
Consequently, 


an > 1, lim an 2 J; 
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We have proved more than stated. It is only the nonva- 
nishing of the limit which is needed. 


Proof. Let’s restate Stirling’s formula this way: 


tim C/V _ 


n= n! 


To prove the result, we need only to show that 


r= lim ap = V2rT. 


n00 

We use the fact that 
; (7h)72?” 
lim -———~ = yT 

n—00 (2n)!/n vr 
which was derived in an earlier section, to evaluate this 
limit r. The function of n appearing here can be rewritten 
in terms of a,, as follows: 

(n!)?2? a? 1 


n 


(n)iyn am V2 


As n becomes infinite, this quotient approaches \/7 on 
the one hand and r?/(rv2) on the other. Hence, r = 
V2, and the proof is complete. 


Corollary. 
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Proof. 


(Qnjle*_ 


lim = lim 
noo (2n)? n>% (oe 


Corollary. 


n+p)! 
Uaa n > œ; p = 1,2,.... 
n! 


Proof. By Stirling’s formula, 


(n+p)! 


ne ™ Inn (1 + (p/n))" P2 
n! eP 


=. 1. 


Each of the three quotients on the right approaches 1. 


Corollary. 
1, 1 
lim —Vn! = =. 
noo n e 


Proof. By Stirling’s formula, 


cio ) Varn = +00 


lim val = lim ( n! n (2rn) m 
(n/e)}} v 2mn 


noo nN n— o0 
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Residues of the Gamma 
Function 


If a complex function is analytic on a region R, it is in- 
finitely differentiable in R.t A complex function may fail 
to be analytic at one or more points through the pres- 
ence of singularities, or through the presence of branch 
cuts. A single-valued function that is analytic in all but 
a discrete subset of its domain, and at those singularities 
goes to infinity like a polynomial (i.e., these exceptional 
points are poles), is called a meromorphic function. In 
this context, the word ” pole” is used to denote a singu- 
larity of a complex function. f has a pole of order n at zo 
if n is the smallest positive integer for which (z— zo)” f(z) 
is analytic at Zo. 


+ An intermediate understanding of complex analysis is neces- 
sary to understand the material in this and the next chapter. 
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Definition 16.1. A function f has a pole at zo if it can 
be represented by a Laurent series centered about zo with 
only finitely many terms of negative exponent, i.e., 


OO 


te) = ` aplz — zo)" 


k=—n 


in some neighborhood of zo, with a-n 4 0, for some n € 
N. The number n is called the order of the pole. A simple 
pole is of pole of order 1. 


Definition 16.2. The constant a_, in the Laurent series 


OO 


f(z) = J alz- z)" 


k=—00 


of f(z) about a point zp is called the residue of f(z). 
The residue of a function f at a point zp is often denoted 
Res(f; zo). 


If f is analytic at zo, its residue is zero, but the converse 
is not always true. 


The residues of a meromorphic function at its poles char- 
acterize a great deal of the structure of a function — 
residues appear in the residue theorem of contour inte- 
gration, which, briefly, is a very powerful theorem used 
to evaluate contour integrals of analytic functions. 
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The residues of a function f(z) may be found without 
explicitly expanding into a Laurent series. If f(z) has a 
pole of order n at zo, then a, = 0 for k < —nanda_, # 0. 


Thus, 


f(z) = X alz- 2)" =Y arne- E 


k=—-n 


(z — 20)" f(z =Y arne a) 
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Differentiating, 
qi lle — 20)" F(] = X kar-nlz = 20)" 
k=0 
= ` kag-n(z — zo)" 7! 
k=1 
=) (k+ 1)ak-n41(2— 2) 
k=0 
a = 
qa [e — 20)"£2)] = J klk + Vanna = 20)" 
k=0 
= X k(k + 1)ar-ny1(2 — 2)** 
k=1 
= Ņ (k + 1)(k + 2)ar-n+2(2 — 20) 
k=0 
Iterating, 
qr-t 


sym e—a)" fe) = SEADE) (+n) eager 


OO 


(n—1)! a_1 b> (k4 -< (k+n—1) ay_i1(z—29)* 


k=1 
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so we have 
qr-t 
im el [(z — 20)" f(2)] = lim (n —1)!a_,+0 
= (n T 1)! a_ 1 


and the residue is 


1 qr} ( 
(n — 1)! dzr-1 


a, = 


a zo)" f (2) 


Z=% 


Proposition 16.1. In the case that the limit lim,_,..(z— 
zo) f(z) exists and has a non-zero value r, the point z = zo 
is a pole of order 1 for the function f and 


Res(f; zo) =r 
This result follows directly from the preceding discussion. 


Now, the Gamma function may be expressed 

T 

i= (z +n) 
z(z+1)\(z +2) (z+n- 1) 

According to the standard definition 


T(z) := f tle dt, 
0 


T(z) is defined only in the right half-plane Re(z) > 0, 
whereas equation 16.1 is defined and meromorphic in the 
half-plane Re(z) > —n where it has poles of order 1 at the 


(16.1) 
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points 0, —1, —2,...,—(n — 1). Equation 16.1 is the ana- 
lytic continuation of T(z) to the half-plane Re(z) > —n. 
Since the positive integer n can be chosen arbitrarily, Eu- 
ler’s Gamma function has been analytically continued to 
the whole complex plane. 


For determining the residues of the Gamma function at 
the poles, rewrite equation 16.1 as 

Tiz+n+1) 
(z+ 1)(2+2)---(z+n) 


r(z) = : 


from which we see 


Pietn ti) 
(z +n)r(z) = z(z+ 1) (z+n-1) 


At the point z = —n, 
r(z+n+1)=T(1)=0!= 
and 
2(z+1)---(¢+n—-1) = (-1)"n! 
which, taken together, imply that 


Cpe 


n! 


Res(T; —n) = 


1/ 


Hankel’s Contour Integral 
Representation 


In this chapter we find integral representations of the 
Gamma function and its reciprocal. Numerical evalua- 
tion of Hankel’s integral is the basis of some of the best 
methods for computing the Gamma function. 


Consider the integral 


T= | (=) edt, 
c 


where Re(z) > 0 and where the contour C starts slightly 
above the real axis at +oo, then runs down to t = 0 
where it winds around counterclockwise in a small circle, 
and then returns to +00 just below the real axis. 
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The ¢-plane is cut from 0 to co. We define the phase such 
that arg(—t) = 0 on the negative real t-axis. On the con- 
tour C we have —7 < arg(t) < a. Therefore, just above 
the positive real t-axis we have arg(—t) = —7, whereas 
just below we have arg(—t) = +7, the angle measured 
counterclockwise. 


Given these definitions, it then follows that 


(ep 1 =e in(z 1) 42 1 


just above the positive real axis, and 


(t) 1 — etin(z 1) 42 1 
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just below the positive real axis. And on the small circle 
enclosing t = 0 we have —t = de’. So then our integral 
becomes 


ô œo 
Te | eiTe) lemt dt + | etim(z-lyzz-le-t dt 
oe) ô 


i / (5e) z—1 eô(cos0+i sin 6) 5619; dé 


T 


> -2isin(rz) f edt for ô— 0. 
0 


The expression above contains the usual representation 
for the Gamma function, so T = —2isin(az)I(z), or 


La re ae dt. (17.1) 


2isin rz 


This is Hankel’s integral representation for the Gamma 


function, valid for all z # 0, +1, +2,.... It has several 
equivalent forms. Considering the fact that e*" = —1, 
we find 


—1 
aa | ea 
2isin Tz Jo 
—1-(-1 z—1 
= ( ) Ta 
C 


2isin rz 


eo ite 
eitz — e-inz c 


1 
OE f the dt, 
emiz — Ije 
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so another form of the contour integral representation is 
given by 


1 
: z= >t ; 
Tea = ae f e “dt 


The trivial substitution s = —t produces yet another 
equivalent form, 


1 
a= azl. s7 te ds. 
2isin Tz Jo 


An application of Euler’s reflection formula, 


T 


r()r(1-—z)= = : 
sin mz 

to Equation 17.1 leads immediately to the contour in- 

tegral representation for the reciprocal of the Gamma 

function, 


1 sin mz 
T(z) T (=g) 
sin —1 ted 
© m Qésin(a(1 — z)) fi EN 
= Sz ot 


So we have for the contour integral representation of the 
reciprocal gamma function, 


l i a t 
voe A et de. 
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An equivalent representation is found by making the sub- 
stitution s = —t, 


1 1 —z sS 
TE =z O ef ds. 


Sometimes we write oe for f c Meaning a path starting 
at infinity on the real axis, encircling zero in the positive 
(counterclockwise) sense and then returning to infinity 
along the real axis, respecting the cut along the positive 
real axis. 


The Weierstrass. Factor 
heorem 


The most important property of the entire rationals is 
expressed in the fundamental theorem of algebra: Every 
non-constant entire rational function has zeros. If we 
have an arbitrary, non-constant entire rational function, 


go(z) = ao + a12 +--+ + amz”, (m > 1, am # 0) 


then it follows from the fundamental theorem of algebra 
that go(z) can be written 


Go(2) = Am(2 = 21)" (z = 22) +++ (z= ZR) 


where 21, 22,..., 2% denote all the distinct zeros of go(z), 
and a1, Q2,...,@% denote their respective orders. 
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For every entire rational function there is a so-called fac- 
tor representation. From this representation, we infer 
that every other entire rational function g(z) which has 
the same zeros to the same respective orders can differ 
from go(z) only in the factor am. It is thus always possible 
to construct an entire rational function (with prescribed 
zeros) that can be represented as a product which dis- 
plays these zeros. 


Suppose the entire function to be constructed is to have 
no zeros at all. Then the constant 1, or the function e7, or 
e”, or more generally, e”) is a solution of the problem, 
if h(z) is a completely arbitrary entire function. 


Theorem A.1. If h(z) denotes an arbitrary entire func- 
tion, then e” is the most general entire function with 
no zeros. 


Proof. We have only to show that if H(z) = ao + az + 
agz7+---+ isa given entire function with no zeros, then we 
can determine another entire function h(z) = bo + biz + 
bəz? +--+» such that e”) = H(z). Now, since H(z) 40, 
we have in particular aj) = H(0) 4 0. Hence, bo can be 
chosen such that e’° = ag; for, e” takes on every value 
except zero. Likewise, HE is everywhere single-valued 
and regular, and is therefore an entire function. 

The same is true of H'(z), so that 


H'(2) 


= Co + C12 + C22 +::: 
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is also an entire function, and this series is everywhere 
convergent. 


The same is true of the series 
z eee 
2 n 
= bo + biz +- + bnz" +e 


which represents an entire function, h(z). If we set e*® = 
H(z), then 


H(z) 
H(z) 


= Co + C12 + C22 as 


and hence H; - H’ — H- H; = 0. Therefore 


H-H -HH df ED) o 
H? de AZ) fe 


and the quotient of the two functions H(z) and H(z) is 
constant. For z = 0 the value of this constant is 1. Thus, 


H(z) = H (2) =e). 


Having thus described the case for which no zeros are 
prescribed, its easy to see the extent to which a general 
entire function is determined by its zeros. If Go(z) and 
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G(z) are two entire functions which coincide in the po- 
sitions and orders of their zeros, then their quotient is 
also an entire function, but one with no zeros. G(z) and 
Go(z) thus differ by at most a multiplicative entire func- 
tion with no zeros. Conversely, the presence of such a 
factor of Go(z) does not alter the positions or orders of 
its zeros. 


Theorem A.2. Let Go(z) be a particular entire function. 
Then, if h(z) denotes an arbitrary entire function, 


G(z) =e") . Go(z) 


is the most general entire function whose zeros coincide 
with those of Go(z) in position and order. 


The question of the possibility and method of construct- 
ing a particular entire function with arbitrarily prescribed 
zeros remains. An entire function has no singularity in 
the finite part of the plane; therefore it can have only 
a finite number of zeros in every finite region. The pre- 
scribed points must not have a finite limit point as a 
consequence. If we make this single restriction, an entire 
function can always be constructed. It can be set up in 
the form of a product which has the positions and orders 
of its zeros. 


Theorem A.3. (Weierstrass factor theorem). Let 
any finite or infinite set of points having no limit point 
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be prescribed, and associate with each of its points a def- 
inite positive integer as order. Then there exists an en- 
tire function which has zeros to the prescribed orders at 
precisely the prescribed points, and is otherwise different 
from zero. It can be represented as a product, from which 
one can read off again the positions and orders of the 
zeros. If Go(z) is one such function, 


G(z) = è . Go(z) 


is the most general function satisfying the conditions of 
the problem, if h(z) denotes an arbitrary entire function. 


The entire function satisfying the conditions of the Weier- 
strass factor theorem will be set up in the form of a prod- 
uct; in general, in the form of an infinite product. As with 
infinite series, we shall assume the reader is familiar with 
the simplest facts in the theory of infinite products. We 
thus present, without proofs, the most important defini- 
tions and theorems for our purposes. 


Definition A.1. The infinite product 
Uy ges ty = [fw (A.1) 
v=1 


in which the factors are arbitrary complex numbers, is 
said to be convergent iff from a certain index on, say for 
all v > m, no factor vanishes, and 


lim (Um+1 *“Um42°°* Un) 
Noo 
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exists and has a finite value distinct from zero. If we call 
this limit Um, then the number 


U =u: U2: Um Um 


which is obviously independent of m, is regarded as the 
value of the infinite product A.1. 


Theorem A.4. A convergent product has the value zero 
iff one of its factors vanishes. 


Theorem A.5. The infinite product A.1 is convergent 
iff having chosen an arbitrary e > 0, an index no can be 
determined such that 


[ent *Unt2°°° Un+r T 1| <E 
for alln > no and allr > 1. 


Since on the basis of this theorem (let r = 1 and n +1 = 
v) it is necessary that lim,...U, = 1, one usually sets 
the factors of the product equal to 1 +c, so that instead 
of dealing with (A.1) one is concerned with products of 


the form 


J[0 +6). (A.2) 


ea 


For these, c, — 0 is a necessary, but insufficient, condi- 
tion for convergence. 
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Definition A.2. The product A.2 is said to be abso- 
lutely convergent if 


OO 


[[c + |e,|) converges. 


v=1 


Theorem A.6. Absolute convergence is a sufficient con- 
dition for ordinary convergence; in other words, the con- 
vergence of [[(1+|c,|) implies convergence of | [(1 + c»). 


On the basis of this theorem it is sufficient for our pur- 
poses to have convergence criteria for absolutely conver- 
gent products. The next two theorems settle the question 
of convergence for such products. 


Theorem A.7. The product | |(1 + Ww), with yw > 0, is 
convergent iff the series X` yy converges. 


Theorem A.8. For [[(1+c) to converge absolutely, it 
is necessary and sufficient that X` c, converge absolutely. 


The following theorem is similar to one on absolutely 
convergent series: 


Theorem A.9. If the order in which the factors of an 
absolutely convergent product occur is changed in a com- 
pletely arbitrary manner, the product remains convergent 
and has the same value. 
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In addition to products with constant factors, we need 
products whose factors are functions of a complex vari- 
able z. We write these products in the form 


OO 


[[G+ Ae). (A.3) 


v= 1 


We designate as the region of convergence of such a prod- 
uct the set M of all those points z which (a) belong to 
the domain of the definition of every f,(z), and for which 
(b) the product A.3 is convergent. Accordingly, the prod- 
uct assumes a certain value for every z of M; thus, the 
product represents in M a certain (single-valued) func- 
tion. For our purposes, it is important to possess useful 
conditions under which such a product, in its region of 
convergence, represents an analytic function. The follow- 
ing theorem suffices: 


Theorem A.10. Let f(z), fo(z),..., fu(z),... be an in- 
finite sequence of functions, and suppose a region G exists 
in which all these functions are regular. Let Y? |fu(z)| 
be uniformly convergent in every closed subregion G’ of G. 
Then the product A.3 is convergent in the entire region 
G, and represents a regular function f(z) in G. More- 
over, by theorem A.1, this function has a zero at those, 
and only those, points of G at which at least one of the 
factors is equal to zero. The order of such a zero is equal 
to the sum of the orders to which these factors vanish 
there. 
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Proof. Let G’ be an arbitrary closed subregion of G. For 
every m > 0, 


XO IA, along with X` |f.(2)| 
v=m+1 v=1 
converges uniformly in G’. By theorem A.8, the product 


OO 


ID a+ fol2)) (A.4) 


v=m4+1 


is absolutely convergent in G’, and represents a certain 
function there. Let us call this function Fm(z). Choose 
m such that 


Pal flH arll (AS) 


for all n > m, all r > 1, and all z in G’ — then F,,(z) is 
regular and nonzero in G’. Indeed, if, for n > m, we set 


n 


lI (1+ f.(z))=P, and P, =0 


v=m4+1 


then we have 
Fm(2) = lim P, 


n— o0 
= lim [(Pmti — Pm) + (Pm+2— Pm+i) +++ (Pa = Pra) 


or 
oo 


F(z) = ` (P, i Ps1) (A.6) 


v=m4+1 
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and Fm(z) is thus represented by an infinite series. Now, 
for n > m, 


|Pal < (1+ |fm+i(z)) A fala) 
< elin OHHO < ef <2 


so the inequality 
|P, R Pl = =A ‘ | fo(z)| < 2| falz) 


is valid for the terms (from the second onward) of the 
series just obtained. Consequently, the new series A.6, 
along with $` |f.(z)|, is uniformly convergent in G’, and 
the function Fm(z) defined by that series is a regular 
function in G’. It is also distinct from zero there. For, by 
(A.5), we have in G’, for n > m, 


RORE 


and hence, for v > m +1, 


1 


ierosina 


so that no factor of Fm can be equal to zero. Since 
f(z) = (1+ filz)) o A finl2)) - Fin(2) 


we have that f(z), together with Fm(z), is regular at ev- 
ery point z of G’, and can vanish at such a point only if 
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one of the factors appearing before Fm(z) vanishes. The 
order of such a zero is then equal to the sum of the orders 
to which these factors vanish there. 


Now let z be an arbitrary point of G. Since z is an interior 
point of G, it is always possible to choose G’ such that z 
also belongs to G’. Hence, the above considerations hold 
for all of G, and the proof is complete. 


We can also make an assertion concerning the derivative 
of f(z). Since the ordinary derivative of a product of 
many factors is difficult to calculate, we introduce the 
logarithmic derivative. We have the following theorem 
concerning this derivative: 


Theorem A.11. Given theorem A.10, we have that 


Me) 5 fie) a 


for every point z of G at which f(z) # 0; i.e., the se- 
ries on the right 1s convergent for every such z and is 
equivalent to the logarithmic derivative of f(z). 


Proof. If z is a particular point of the type mentioned 
above, and if the subregion G’ is chosen so as to contain 
z, then 


Pz) _ O inl), Fal) a) 


O tee). Pe) EL): 
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Since the series A.6 converges uniformly in G’, 


/ = I1 pl ee 1 
Falz) = 2R P24) tim P,, 


according to theorem A.6. Here P’ denotes the derivative 
of Pa. Since F;,(z) and all P, for n > m are not zero, 


OT: E ee AC) 
Falz) o nao Pa no 1+ fm+i(2) | 1+ falz) 
wa 
2s 1+ fol) 


which, with (A.8), proves the assertion. 


Theorem A.12. The series (A.7) converges absolutely 
and uniformly in every closed subregion G" of G contain- 
ing no zero of f(z), and hence may be repeatedly differ- 
entiated there any number of times term by term. 


Proof. Since none of the factors (1 + f,(z)) can vanish 
in G”, the absolute value of each remains greater than 
a positive bound, y, say. Since this is certainly greater 
than 1/2 for all v > m, a positive number y exists, such 
that yy > y > 0 for all v. Then, for all v and all z in G”, 


fi) 
1+ folz) 
From the proof of theorem A.6 it follows that S>|f/(z)| 


converges uniformly in G”. By the last inequality, this is 
also true then of the series A.7. 
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Having become familiar with infinite products, it is straight- 
forward to prove the Weierstrass factor theorem. 


Proof. If only a finite number of points 21, 22,..., Zk hav- 
ing the respective orders a 1,Q2,...,Q@,% are prescribed, 
then the product 


(z — 21) (z — 22) +++ (z — ZR) (A.9) 


is already a solution of the problem, so that this case is 
settled immediately. If, however, an infinite number of 
points are prescribed as zeros, because the corresponding 
product would be meaningless generally. This would still 
be the case if, with regard to the infinite products we 
are dealing with here, we were to replace (A.9) by the 
product 


(1-2 (1-2 (1-2) (A.10) 


which serves the same purpose. We therefore modify our 
approach — and therein lies the utility of Weierstrass’s 
method. 


The set of prescribed points is enumerable, since every 
finite region can contain only a finite number of them. 
They can therefore be arranged in a sequence. The way 
in which the points are numbered is not important. How- 
ever, if the origin, with the order ag, is contained among 
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them, we shall call this point zo and, leaving it aside for 
the present, arrange the remaining points in an arbitrary, 
but then fixed, sequence: 21, 22,...,Zy,.... Let the cor- 
responding orders be ay, Q2,...,Qy,.... The z are all 
different from zero; and since they have no finite limit 
point, 

Zy > 00, [æl 4 +. 


Consequently, it is possible to assign a sequence of posi- 
tive integers k1, ko,...,ky,... such that 


oo ky 

z 

So (=) (A.11) 
v=1 

is absolutely convergent for every z. In fact, it suffices, 

e.g., to take k, = v +a. For, no matter what fixed value 


z may have, since z, — oo, we have for all sufficiently 
large v 


and hence 


2 1 lace 1 v 
"RD 2 


and the absolute convergence of the series is thus assured. 


129 


Let the numbers k, be chosen subject to this condition, 
but otherwise arbitrarily, and keep them fixed. Then we 
shall prove that the product 


Gol2) = TT (0 = 2) opl +: (y eee et (yyy 


represents an entire function with the required proper- 
ties. (Here the factor z™ appearing before the product 
symbol is to be suppressed in case the origin is not one 
of the prescribed zeros. Likewise, if one of the numbers 
ky, is equal to unity, the corresponding exponential factor 
does not appear.) 


The proof of this assertion is now very simple. To be able 
to apply our theorems on products, we set the factors of 
our infinite product equal to 1+ f,(z). According to 
theorem A.10, we must then merely prove that 


EROSEN (: - z) 


(A.12) 


ky—1 Qy 
z 1 z 1 
-ex | Meso | = 
p A ky -—1\ 2, 


converges uniformly in every bounded region. For then 
the entire plane can be taken as the region G of theorem 
A.10, according to which the infinite product, and con- 
sequently also Go(z), is an entire function. On account 
of the form of the factors of Go(z), the second part of 
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theorem A.10 at once yields that Go(z) also possesses 
the required properties. The uniform convergence of the 
series (A.12) in the circle about the origin with radius R 
(R > 0 arbitrary, but fixed) is established as follows: 


Since the series (A.11) also converges for z = R, and 
since zy —> oo, m can be chosen so large that 


(A.13) 


for all v > m. Let us for the moment replace z/z, by u, 
ky by k, and a, by a. Then, for v > m, the vth term of 
the series (A.12) has the form 


[o-oo bs Sate = Ti 


with |u| < 4 and alu|* < £. Now for |u| < 1 we can set 


ae. ae 
EEE S rare a e wen en 
u exp { u 5 3 \ 


so that this vth term is further equal to 
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and hence 


7 |u|” | juj +t | =i 
RRNA E 


< exp fauta + Jul + ul? +- } sige ia] 


because |u| < E, Further, since e” — 1 < xe” for x > 0, 
the vth term is less than or equal to 


Saluer e gauk 


the exponent of e being smaller than one, according to 
(A.13). Hence, for all sufficiently large v and all |z| < R 
we have 


ky ky 


Ifo(2)| < 6a|— 


v 


R 
< 6a, oan 
ey 


But these are positive numbers whose sum converges (be- 
cause of the manner in which the k, were chosen). There- 
fore, by Weierstrass’s M-test, |f,(z)| is uniformly con- 
vergent in the circle with radius R about the origin as 
center, so the proof of the Weierstrass factor theorem is 
complete. 


The product is simplest if the prescribed zeros and orders 
are such, that the series }>a@,,/z,, and consequently, for 
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every z, the series $` a,(z/zv), converges absolutely for 
our sequence 21, 22,.... For then it is possible to take all 
k, = 1, and the desired function is obtained simply in 
the form 


If, for example, the points 0,1,4,9,...,v?,... are to be 
zeros of order unity, then we have that 


G(z) = eM). z. I] | — 3 
v=1 


with h(z) an arbitrary entire function, is the most general 
solution of the problem. If the points 1,8,...,v°,... are 
to be zeros of respective orders 1,2,...,v,..., then 


h(z = Z 
Gz) = t. (2-3) 
v=1 


is the most general solution. 


Now we present an application of the factor theorem 
which is of importance to an earlier topic in the text, 
namely, in the derivation of Euler’s reflection formula 
where we used an identity known as the product repre- 
sentation of the sine function. 
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Consider constructing an entire function which has zeros, 
of order unity, at precisely all the real lattice points (i.e., 
at 0,+1,+2,...). We number these points so that zp) = 
Os 29 = +1, eo = —1,...,Z%w-1 = Veoh = —v,..., with 
(v = 1,2,...). The series 


v=1 


is absolutely convergent for every z, and we can therefore 
take all k, = 2. 


Then we have 


G(z) = eb. z. II (2 — Jem 


ome II (1 = Jer 
v=1 g 
=¢.z || (1-5) 


as the most general solution of the problem. 


Since the function sin mz is evidently also a solution of 
the problem, it must be contained in the expression just 
found. That is, there exists a certain entire function, 
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which we shall call ho(z), such that 


sin mz = eh) . z. I (: — 5) (A.14) 


If we can succeed in obtaining this function ho(z), we 
shall have the factor representation of sin mz. 


The function ho(z) certainly cannot be ascertained from 
a knowledge of the zeros alone. On the contrary, for its 
determination we must make use of further properties of 
the function sin 7z; e.g., its power series expansion, its 
periodicity properties, the conformal map effected by it, 
its behavior at infinity, and so on. Now to determine 


ho(z). 


First, to show that hj(z) is a constant. According theo- 
rem A.11, it follows from (A.14) that 


1 G/ 1 1 
nm cot mz = holz) + A ( + j (A.15) 
v=1 


Z—-v z+v 


According to theorem A.12, this expression may be dif- 
ferentiated repeatedly term by term. Thus, 


T 1 
“eee 8 9-3-5 7 y) 
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or, more succinctly, 


ne ja y 1 T’ 
DAS (z= v)? sin?rz 


v=— o0 


This relation holds in every closed region which contains 
no real lattice points. If we replace z by z + 1 in the 
rightmost member, it is not altered because sin?7z has 
the period +1, and 


= 1 = 1 ad 
2 (2+1-vp 7È: (2—(— Dp os (z= nw)? 


Hence, hj(z) is an entire function with the period +1. 
In order to show that h§(z) is a constant, it is sufficient 
to show that |hj(z)| cannot become arbitrarily large. On 
account of the periodicity of hg(z) which we just estab- 
lished, it is sufficient, for this purpose, to show that a 
constant K exists such that |hg(z)| < K for all z = z+ iy 
for which 0 < x < 1 and |y| > 1. 


Now for these z, 


+00 1 +00 1 [0:0] 1 
ey ON 
2 (z= v)? <>, (x= v) +y Pee 


and, since sin mz = (1/27)(e’"* — e~*””), 


2 


1 E An? An? 


= < 
e27y + e277 — 2cos2rz ` ell —2 


sin? rz 
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for those z. Consequently, 


= 1 Ar? 
ooe ! e27lyl — 2 


there, and this expression certainly remains less than a 
fixed bound for all |y| > 1. Hence, 


ho(2) 


ho(z) = constant = c. 


According to the inequality just obtained, |hg(z)| is arbi- 
trarily small if |y| is sufficiently large; hence c” must be 
equal to zero. Therefore 


holz) =0, ho(z) = constant = c 
and hence by A.15 
2z 


1 
Tcotmz=c 4 X 5 z: 
z zł — v 


v=1 


Now if we substitute —z for z in this equality, we see that 
cd = —c’, and hence c = 0. Then ho(z) and e”?™ are also 


constant. Therefore 


If we divide through by z and allow z to approach zero, 
we obtain 7 = c. We thus have, 


a 2 
mesae ll I-35 
valid for all z, which is the product representation of the 
sine function we set out to find. 


The Mittag-Leffler Theorem 


For every fractional rational function there is a so-called 
decomposition into partial fractions, in which the poles 
and corresponding principal parts are apparent. Thus, let 
fo(z) be the given rational function, and let 21, z2,..., Zk 
be its poles with the corresponding principal parts 


(v) (v) (v) 
a a a 
hy(z) = — + — 5 te 2 B.1 
(2) Z—By (2-24)? (z — Zy)% ey 
where v = 1,2,...,k. Then we can set 


folz) = go(z) + hi(z) + halz) + +++ + Ag (z) (B.2) 


where go(z) is a suitable entire rational function. From 
this decomposition into partial fractions, we infer that 
every other rational function f(z) having the same poles 
with the same respective principal parts can differ from 
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fo(z) in the term go(z) alone. Furthermore, one can ar- 
bitrarily assign these poles and their principal parts. In 
other words, it is always possible to construct a rational 
function whose poles and their principal parts are pre- 
scribed. This function can be represented as a partial- 
fractions decomposition which displays these poles and 
their principal parts. The most general function of this 
kind is obtained from a particular one by adding to it an 
arbitrary entire rational function. 


These fundamentals concerning rational functions can be 
carried over to the more general class of meromorphic 
functions. 


Definition B.1. A single-valued function shall — without 
regard to its behavior at infinity — be called meromorphic, 
if it has no singularities other than at poles in the entire 
plane. 


Theorem B.1. A meromorphic function has in every 
finite region at most a finite number of poles. 


For otherwise there would exist a finite limit point of 
poles, and this point would be singular, but certainly not 
a pole. 


According to this, the rational functions are special cases 
of meromorphic functions, and the entire functions must 
also be regarded as such. 
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The function 1/sin z is meromorphic because in the fi- 
nite part of the plane it has a singularity, namely a pole 
of order unity, only wherever sin z has a zero. We see, 
likewise, that cot z = cos z/ sin z and tan z are meromor- 
phic functions. More generally, if G(z) denotes any entire 
function, its reciprocal, 1/G(z), is a meromorphic func- 
tion. For, 1/G(z) has poles (but otherwise no singulari- 
ties) at those, and only those, points at which G(z) has 
zeros; and the orders of both are the same. If Gj(z) is an 
entire function which has no zeros in common with G(z), 
we see that Gi(z)/G(z) is a meromorphic function whose 
poles coincide in position and order (although, in general, 
not in their principal parts!) with those of 1/G(z). 


We inquire whether, and how, one can construct a mero- 
morphic function if its poles and the corresponding prin- 
cipal parts are prescribed, and to what extent a mero- 
morphic function is determined by these conditions. 


If Mo(z) and M(z) are two meromorphic functions which 
coincide in their poles and the corresponding principal 
parts, then their difference, M(z)—Mo(z), is evidently an 
entire function. Consequently, they differ by at most an 
additive entire function (a meromorphic function with no 
poles). Conversely, since the addition of such a function 
to Mo(z) does not alter its poles or the corresponding 
principal parts, we are able to say: 


Theorem B.2. Let Mo(z) be a particular meromorphic 
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function. Then, if G(z) denotes an arbitrary entire func- 
tion, 

M(z) = Mo(z) + G(z) 
is the most general meromorphic function which coin- 
cides with Mo(z) in its poles and the corresponding prin- 
cipal parts. 


There remains only the possibility and method of con- 
structing a particular meromorphic function with arbi- 
trarily prescribed poles. 


According to theorem B.1, the set of assigned poles can- 
not have a finite limit point. If this is excluded, however, 
then the problem posed can be solved without any fur- 
ther restriction. 


Theorem B.3. (Mittag-Leffler partial-fractions the- 
orem). Let any finite or infinite set of points having no 
finite limit point be prescribed, and associate with each 
of its points a principal part, i.e., a rational function of 
the special form (B.1). Then there exists a meromor- 
phic function which has poles with the prescribed prin- 
cipal parts at precisely the prescribed points, and is oth- 
erwise regular. It can be represented in the form of a 
partial-fractions decomposition from which one can read 
off again the poles along with their principal parts. Fur- 
ther, by theorem B.2, if Mo(z) is one such function, 


M(z) = Mo(z) + G(z) 


is the most general function satisfying the conditions of 
the problem, if G(z) denotes an arbitrary entire function. 
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If we let M(z) be an arbitrarily given meromorphic func- 
tion, the set of its poles has no finite limit point. Hence, 
according to Mittag-Leffler’s theorem, another meromor- 
phic function, M,(z), having the same poles and prin- 
cipal parts as M(z), can be constructed in the form of 
a partial-fractions decomposition displaying these. Then 
by theorem B.2, 


M(z) = Molz) + Go(z) 


where Go(z) denotes a suitable entire function. We have 
thus actually obtained a decomposition of the given mero- 
morphic function M (z) into partial fractions, from which 
its poles and the corresponding principal parts can be 
read off. 


Proof. If the function to be constructed is to have no 
poles at all, then every entire function is a solution of 
the problem. If it is to have the finitely many poles 
21, 22,---, 2 With the respective principal parts h,(z), 
ho(z), ..., Ae(z), then 


Molz) = hi(z) + haz) +--+ + Ae (2) 


is a solution. If, however, an infinite number of poles are 
prescribed, we cannot attain our goal simply because the 
analogous series, being infinite, would generally diverge. 
Nevertheless, we can produce the convergence by means 
of a suitable modification of the terms of the series. 
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If the origin is a prescribed pole, we denote it by z and 
leave it aside for now. Let ho(z), hi(z),...,Rv(z),... be 
the principal parts corresponding to the points Zo, 21, ..., 


Zy, --.; h,(z) is understood to be an expression of the 
type appearing in (B.1). Each of these functions h,(z), 
v = 1,2,3,..., is regular in a neighborhood of the origin. 


Its power-series expansion 
hy(z) =a + Mz 4a? +... (v =1,2,...) 


for this neighborhood converges for all |z| < |z,|; hence, 
it is uniformly convergent for all |z| < ł|z„|. Conse- 
quently (for every v = 1,2,3,...) an integer n, can be 
determined such that the remainder of the power series 
after the n,th term remains, in absolute value, less than 
any preassigned positive number, e.g., 1/2”. Denote the 
sum of the first n, terms of the series by g,(z). Thus, 
gv(z) is an entire rational function of degree n,: 


galz) = aÙ | az peg (= 1D Senne) 
and for all |z| < 3|z,| we have 
1 
Inu (2) = a2) < go. 


Then 


Mo( ) + Doth )— ghz z) 
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is a meromorphic function satisfying the conditions of the 
theorem. (If the origin is not assigned as a pole, the term 
ho(z) must be omitted. 


To prove this, we must merely show that the right-hand 
side defines an analytic function having in every finite 
domain, e.g., a circle with radius R about the origin as 
center, exactly the prescribed singularities and no others. 


Now, |z,| + +oo. Therefore it is possible to choose m so 
large, that |z,| > 2R, and hence R < $|2,|, for all v > m. 
Then, for all |z| < R and all v > m, 


1 1 
|z| < zll and consequently  |he(z) — gu(z)| < —. 


2v 


Hence, for all |z| < R, the series 


oO 


` [ha (z) = 9v(Z)] 


v=m4+1 


is absolutely and uniformly convergent. Since its terms 
are regular for |z| < R (because the poles of the h,(z) 
with v > m lie outside the circle |z| = R), it defines there 
a regular function which we shall denote by F;,,(z). Then 
evidently 
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is also an analytic function which is regular in the circle 
with radius R about the origin as center, with the ex- 
ception of those points z, in this circle which are poles 
with principal parts h,(z). The same is valid for every 
finite region, because R was completely arbitrary and 
hence, Mo(z) is a meromorphic function with the neces- 
sary properties. 


From the proof it follows that it is sufficient to take the 
degree n, of the polynomial g,(z) (the sum of the first 
Ny, terms of the power series for h,(z)) so large that hav- 
ing chosen an arbitrary R > 0, the terms |h,(z) — g.(z)| 
for all |z| < R finally (i-e., for all sufficiently large v) re- 
main less than the terms of a convergent series of positive 
terms. 


The convergence terms g,(z) are not always necessary. 
Then, of course, the function to be constructed is espe- 
cially simple. If, e.g., the points 0,1,4,...,v?,... are to 
be poles of order unity with respective principal parts 
1/(z — v?), then 


Ten | Za 
Mola) Se) a ae 
v=1 v=0 


is a solution. For, let R > 0 be chosen arbitrarily, and 
m > V2R. Then the series from v = m + 1 on is evi- 
dently uniformly convergent in |z| < R, which proves the 
assertion. 
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Consider the case of cot tz. The real lattice points are to 
be poles of order unity with the residue +1, and hence, 
with the principal parts 


1 
hy(z) = gue ’ (Zo = 0, 2%—1 = V, Zw = —v). 


For v = 1,2,3,..., 


because then for all sufficiently large v (namely, for all 
v > 4R) and all |z| < R, 


R P 2R 
lol — R) [z]? 


|hy(z) — J(2)| < [z |( 


so that the |h,(z) — g,(z)| finally remain less than the 
terms of an obviously convergent series of positive terms. 
Consequently, according to the concluding remark of the 
preceding paragraph, if G(z) is an arbitrary entire func- 


tion, 
M(z) =G(z) +5 S a 
=G) +7 S(j i) 
SGi : xi ! = 


is the most general function of the kind required. The 
function cot mz also has poles of order unity at the points 


0,+1,+2,.... If n is one of them, the residue at this 
point is 
lim Za Mcosme im (@= mI" t- J _ 1 
zn sin mZ zn (—1)"1(z = n) Eey T 


which can be read off from the series expansion for a 
neighborhood of the point z = n. Hence, the function 
mcot mz is contained among the functions M(z) which 
we constructed. 


The still undetermined entire function G(z), which there 
was called ho(z), cannot be ascertained solely from the 
nature and position of the poles. We should, as before, 
have to make use of special properties of the function. 
However, in determining the product sin mz, we have al- 
ready discovered that we have to set ho(z), that is, G(z), 
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equal to zero. Therefore 


it, 1 1 
tmz =-4 ` | 
T COU TZ F BS zr 


v=l 


which is the partial-fractions decomposition of the cotan- 
gent function. 


